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Summary. We review angle resolved photoemission spectroscopy (ARPES) re- 
sults on the high Tc superconductors, focusing primarily on results obtained on 
the quasi- two dimensional cuprate Bi2Sr2CaCu2 08 and its single layer counterpart 
Bi2Sr2Cu06. The topics treated include the basics of photoemission and method- 
ologies for analyzing spectra, normal state electronic structure including the Fermi 
surface, the superconducting energy gap, the normal state pseudogap, and the elec- 
tron self-energy as determined from photoemission lineshapes. 

1 Introduction 

Angle resolved photoemission spectroscopy (ARPES) has played a major role in 
the elucidation of the electronic excitations in the high temperature cuprate super- 
conductors. Several reasons have contributed to this development. First, the great 
improvement in experimental resolution, both in energy and momentum, aided by 
the large energy scales present in the cuprates, allows one to see features on the 
scale of the superconducting gap. More recently the resolution has improved to such 
an extent, that now features in traditional superconductors like Nb and Pb, with 
energy scales of a meV, can be observed by ARPES [Q]. 

Second, most studies have focused on Bi2Sr2CaC'u20s (Bi2212) and its sin- 
gle layer counterpart, Bi2Sr2CuO(i (Bi2201). These materials are characterized by 
weakly coupled BiO layers, with the longest interplanar separation in the cuprates. 
This results in a natural cleavage plane, with minimal charge transfer. This is cru- 
cial for ARPES, since it is a surface sensitive technique; for the photon energies 
typically used, the escape depth of the outgoing electrons is only of the order of 
~ lOA. For this reason, we elect to concentrate on these materials in the current 
article, since the data are known to be reproducible among the various groups. 

The third reason is the quasi-two dimensionality of the electronic structure of 
the cuprates, which permits one to unambiguously determine the momentum of the 
initial state from the measured final state momentum, since the component parallel 
to the surface is conserved in photoemission. Moreover, in two dimensions, ARPES 
directly probes the the single particle spectral function, and therefore offers a com- 
plete picture of the many body interactions inherent in these strongly correlated 
systems. 

There are other reasons to point out as well. For example, the very incoherent 
nature of the excitations near the Fermi energy in the normal state does not allow 
the application of traditional techniques, such as the de Haas van Alphen effect. 
The electrons simply do not live long enough to complete a cyclotron orbit. Another 
important reason why ARPES has played a major role is the highly anisotropic 
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nature of the electronic excitations in the cuprates, which means that momentum 
resolved probes are desirable when attempting to understand these materials. 

It is satisfying to compare the earlier ARPES papers to more recent ones, and 
observe the remarkable progress the field has undergone. But this satisfaction must 
be tempered by the knowledge that the current literature is still full of unresolved 
issues. Undoubtedly, some of what is described here will be viewed in a new light 
in the future, and we will try to point out the contentious points, but the reader 
must keep in mind that the writers, by necessity, approach this task with their own 
set of biases. For a complementary view, the reader is refered to the recent review 
by Damascelli, Shen and Hussain An earlier review of ARPES studies of high 
Tc superconductors can be found in the book by Lynch and Olson j^. 

2 Basics of Angle-Resolved Photoemission 

The basics of angle-resolved photoemission have been described in detail in the 
literature[^. We will limit ourselves to a brief review of some salient points, empha- 
sizing those aspects of the technique which will be useful in understanding ARPES 
studies of high Tc superconductors. We start by looking at a simple independent 
particle picture, and subsequently include the effects of strong interactions. 

The simplest model of ARPES is the three step model [^], which separates 
the process into photon absorption, electron transport through the sample, and 
emission through the surface. In the first step, the incident photon with energy hv 
is absorbed by an electron in an occupied initial state, causing it to be promoted 
to an unoccupied final state, as shown in Fig. There is conservation of energy, 
such that 



hu^ BE + 'P + Ei.in (1) 

where <f is the work function, and BE and -Ekin are the binding and kinetic energies 
of the electron, respectively. We will not discuss the subsequent two steps of the 
three step model, photoelectron transport to the surface and transmission through 
the surface into the vacuum, as they only affect the number of emitted electrons, 
and thus the absolute intensity 

The kinetic energy of the electrons is measured by an electron energy analyzer. 
If the number of emitted electrons is then plotted as a function of their kinetic 
energy, as shown in Fig. |l| peaks are found whenever an allowed transition takes 
place. Eq. (1) then yields the kinetic energy of the electron if the work function is 
known. Measuring <f accurately is not a trivial task, but fortunately, in the case of 
a metallic sample one does not need to know By placing the sample in electrical 
contact with a good metal (e.g., polycrystalline gold) one can measure the binding 
energy in the sample with respect to the chemical potential (Fermi level Ef) of Au. 
The photoemission signal from Au will simply be a Fermi function convolved with 
the experimental energy resolution, and from the mid point of its leading edge one 
estimates Ef. 

The existence of the sample surface breaks (discrete) translational invariance 
in the direction normal to the surface. However, two-dimensional translational in- 
variance in the the directions parallel to the surface is still preserved, and thus kn , 



Photoemission in the High Superconductors 



5 





fvs f> -: (1 1 (t 

l-.iiLT*;} ) 



Fig. 1. a: Independent particle approximation view of ARPES; To the right is the 
resulting energy distribution curve (EDC); b: Experimental EDCs along a path in 
momentum space in Bi2212; c: Intensity map from these data. 



the component of the electron momentum parallel to the surface, is conserved in 
the emission process. This allows us to obtain the in-plane momentum of the initial 
state by identifying it with the parallel momentum of the emitted electron 



kn — \ z sm 6 cos c 



(2) 



of the outgoing electron emitted along the {6, 0) direction with a kinetic energy 
i?kin, as shown in Fig. |^. 

The momentum perpendicular to the sample surface is not conserved, and 
thus knowledge of the final state k( does not permit one to say anything useful 
about the initial state kl, except in the case where ky — 0, in which case k^ — k^ . 
We should emphasize that a given outgoing electron corresponds to a fixed, but 
apriori unknown, kl. In materials with three-dimensional electronic dispersion it is, 
of course, essential to fully characterize the initial state, and many techniques have 
been developed to estimate see Ref. Aside from making ky = as described 
above, one can vary the incident photon energy hv thereby changing the kinetic 
energy of the outgoing electron. For fixed ky this amounts to changing the initial 
state kz. (One also has to take into account the changes in intensity brought about 
by the photon energy dependence of the matrix elements). 

We will not be much concerned about methods of fc^ determination here, since 
the high Tc cuprates are quasi-two-dimensional (2D) materials with, in many cases 
(e.g., Bi2212), no observable fcz -dispersion. In fact, this makes ARPES data from 
2D materials much easier to interpret, and is one of the reasons for the great success 
of the ARPES technique for the cuprates. In the remainder of this article, we will 
use the symbol k to simply denote the two-dimensional momentum parallel to the 
sample surface, unless explicitly stated otherwise. 
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In the independent particle approximation, the ARPES intensity as a function 
of momentum k (in the 2D Brillouin zone) and energy uj (measured with respect 
to the chemical potential) is given by Fermi's Golden Rule as 



where ijji and ipf are the initial and final states, p is the momentum operator, and 
A the vector potential of the incident photon ||6|. Here /(w) = l/[exp(w/T) + 1] is 
the Fermi function at a temperature T in units where R = fcs = 1. It ensures the 
physically obvious constraint that photoemission only probes occupied electronic 
states. 

For noninteracting electrons, the emission at a given k is at a sharp energy 
corresponding to the initial state dispersion. As we will discuss below, the effect 
of interactions (either electron-electron or electron-phonon) is to replace the delta 
function by the one-particle spectral function, which is a non-trivial function of 
both energy and momentum. Although now the peaks will in general shift in en- 
ergy and acquire a width, the overall intensity is still governed by the same matrix 
element as in the non-interacting case. Thus the lessons learned from studying gen- 
eral properties of the matrix elements are equally applicable in the fully interacting 
case. 

2.1 Matrix Elements and Selection Rules 

The dipole matrix element Mfi = (i/'/l A ■ p j?/;;) is in general a function of the the 
momentum k, of the incident photon energy hv, and of the polarization A of the 
incident light. Here we discuss dipole selection rules which arise from very general 
symmetry constraints imposed on Mfi, and which are very useful in interpreting 
ARPES data. Our approach is a simplified version of Hermanson's analysis 

Let the photon beam be incident along a plane of mirror symmetry of the sample 
{M), with detector placed in the same mirror plane; see Fig. ga). The final state 
V'/ must then be even with respect to reflection in M, because if it were odd the 
wave function would vanish at the detector. (For this discussion it is simpler to 
imagine changing the photon polarization, keeping the detector fixed. In the actual 
experiment, however, it is the polarization which is fixed and the detector is moved 
when checking dipole selection rules.) 

The dipole transition is allowed if the entire matrix element has an overall even 
symmetry. Thus two possibilities arise. (1) For an initial state ipi which is even with 
respect to M, the light polarization A must also be even, i.e. parallel to M. (2) 
For an initial state odd with respect to M, A must also be odd, i.e. perpendicular 
to M. This can be summarized as: 



Consider hybridized Cu3d — 02p initial states, as shown in Fig. y(b), which 
have a d^2_y2 symmetry about a Cu site. These states are even under reflection 



I{k,u;) cx KV/I A ■ p 1^,)!" f{ij)S{cj ~ et) 



(3) 




(4) 
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Fig. 2. (a) Arrangement of the photon beam and detector in order to make use 
of the photoemission selection rules, (b) Symmetry of the Cu3d^2_y2 orbitals hy- 
bridized with the 02p orbitals. (c) EDCs showing the symmetry of the orbitals in 
(b) obtained at hv =22 eV. 



in (0,0) — (7r,0) (i.e. the plane defined by this symmetry axis and the z-axis) and 
odd with respect to (0, 0) — (tt, tt). Therefore, measurement along the (0, 0) — (tt, 0) 
direction will be dip ole- allowed (forbidden) if the polarization vector A is parallel 
(perpendicular) to this axis. Fig. |2|(c) shows that, consistent with an initial state 
which is even about (0,0) — (7r,0), the signal is maximized when A lies in the 
mirror plane and minimized when A is perpendicular to this plane. (The reasons 
for small non-zero intensity in the dipole forbidden geometry are the small, but 
finite, k-window of the experiment and the possibility of a small misaUgnment 
of the sample). Similarly, we have checked experimentally (for Bi2212 in the Y- 
quadrant where there are no superlattice complications) that the initial state is 
consistent with odd symmetry about (0,0) — {tt,tt) (see also Ref. Q). 

While the dipole matrix elements are strongly photon energy dependent, the 
selection rules are, of course, independent of photon energy. This has been checked 
by measurements at 22 eV and 34 eV.j^ All of these results are consistent with 
the fact that we are probing Cu3d — 02p initial one-electron states with d^2_y2 
symmetry. We will see below that these selection rules are extremely useful in dis- 
entangling the main Cu02 "band" from its umklapp images due to the superlattice 
in Bi2212. 



2.2 The One-Step Model 

The three-step model of photoemission gives a very useful zeroth order description 
of the photoemission process, but it needs to be put on a firmer footing, both as 
regards the calculation of photoemission intensities from ab-initio electronic struc- 
ture calculations, and also for developing a deeper understanding of how many-body 
effects influence ARPES spectra. We now briefly review the so called "one-step" 
model of photoemission with these goals in mind. 
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Ever since Hertz and Hallwachs discovered photoemission p^] , it is known that 
the photo-electron current at the detector is proportional to the incident photon 
flux, i.e., to the square of the vector potential. Thus photoemission measures a 
nonlinear response function, and the relevant correlation function is a three cur- 
rent correlator, as first emphasized by Schaich and Ashcroft It is instructive 
to briefly review their argument. As in standard response function calculations, 
let's look at an expansion of the current at the detector (the response) in powers 
of the vector potential of incident photons (the perturbation). Let R be the loca- 
tion of the detector in vacuum and r denote points inside the sample. The zeroth 
order contribution {0|ja(R, t)jO) vanishes as there are no currents flowing in the 
ground state |0) of the unperturbed system. The linear response also vanishes, with 
{0\ja(R;t)j3{r,t')\0) = and {0\ja{r,t')jp{'R.,t)\0} = 0, since there are no particles 
at the detector, in absence of the electromagnetic field, and j,3(R, t)]0) = 0. Thus 
the leading term which survives is 



(j^(R,t))oc / drdt'dr"dt"Ac{r,t')Af,{r",t") 

{0\Mr',t')j^{Il,t)j,3{v",t"m (5) 



where only current operators inside the sample act on the unperturbed ground state 
on either side and the current at the detector is sandwiched in between. 

The three current correlation function can be represented by the triangle di- 
agram shown in Fig. ^A) where the line between the two external photon 
vertices is the Greens function of the "initial state" or "photo-hole" and the two 
lines connecting the photon vertex to the current at the detector represent the 
"photo-electron" which is emitted from the solid. There is a large literature on 
the evaluation of the (bare) triangle diagram incorporating the results of ab-initio 
electronic structure calculations of a semi-infinite solid, including the effects of re- 
alistic surface termination and multiple scattering effects in the photo-electron final 
states. Detailed calculations using this formalism were first carried out by Pendry 

for a system with one atom per unit cell, and later generalized to more complex 
crystals @,0,[6). 

This approach is very reliable for calculating photoemission intensities, and 
gives important information about "matrix element effects", i.e., the dependence 
of the ARPES intensity on k, and on the incident photon energy and polarization. 
(It does not, however, shed light on the many-body aspects of the lineshape). Such 
ab-initio methods have been extensively used by Bansil, Lindroos and coworkers 
p^ ] for the high Tc cuprates. 

As an example of the usefulness of ab-initio methods, we show the comparison 
between the observed and calculated ARPES intensities for YBCO in Fig. ^. Since 
the calculated intensity is a sensitive function of the termination plane, such a com- 
parison suggests that the crystal breaks at the chains, a conclusion later confirmed 
by STM measurements In fact, it was the complexity of the termination in 
YBCO, plus the fact that none of the possible cleavage planes of YBCO are charge 
neutral, and thus involve significant charge transfer, that convinced us to focus 
primarily on the Bi-based cuprates (Bi2212 and Bi2201). These materials have two 
adjacent BiO planes which are van der Waals bonded, and thus cleavage leads to 
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neutral surfaces which are not plagued by uncontrolled surface effects. The highly 
two-dimensional Bi-based cuprates are thus ideal from the point of view of ARPES 
studies. 

2.3 Single-Particle Spectral Function 

Although the one step model gives a reasonable interpretation of the overall in- 
tensity in the photoemission process, much of what we really want to learn from 
the experiment relates to the spectral lineshape, which, as we will show, is strongly 
influenced by correlations. We must then ask ourselves the important question of 




Fig. 4. Comparison of the photoemission intensity of states along the zone diag- 
onal in YBCO with (a) even polarization and (b) odd polarization, showing that 
this state is odd with respect to reflection about the diagonal to the Cu-O bond 
direction, (c) Comparison of the photoemission intensity to the observed one as a 
function of polarization. The detailed polarization comparison gives confidence to 
the interpretation of the observed peak as arising from the Cu-O plane states. 

how the initial state lineshape enters the ARPES intensity and to what extent 
this is revealed in the observed lineshape. Although a fully rigorous justification of 
a simple interpretation of ARPES spectra is not available at the present time, a 
reasonable case can be made for analyzing the data in terms of the single particle 
spectral function of the initial state. 

We begin by considering the various many body renormalizations of the bare 
triangle diagram which are shown in Fig. ^ These renormalizations can arise, in 
principle, from either electron-phonon or electron-electron interactions. These self 
energy effects and vertex corrections are easy to draw, but impossible to evaluate in 
any controlled calculation. Nevertheless, they are useful in obtaining a qualitative 
understanding of the the various processes and estimating their relative importance. 
Diagram (B) represents the self-energy corrections to the occupied initial state that 
we are actually interested in studying; (C) and (D) represent final state Une-width 
broadening and inelastic scattering; (E) is a vertex correction that describes the 
interaction of the escaping photo-electron with the photo-hole in the solid; (F) is a 
vertex correction which combines features of (D) and (E). (An additional issue in a 
quantitative theory of photoemission is related to the modification of the external 
vector potential inside the medium, i.e., renormalizations of the photon line. These 
are considered in detail in Ref. jl^). 

If the sudden approximation is valid, we can neglect the vertex corrections: 
the outgoing photo-electron is moving so fast that it has no time to interact with 
the photo-hole. Let us make simple time scale estimates for the cuprates with 15 
- 30 eV (ultraviolet) incident photons. The time t spent by the escaping photo- 
electron in the vicinity of the photo-hole is the time available for their interaction. 
A photoelectron with a kinetic energy of (say) 20 eV has a velocity v = 3 x 10** 
cm/s. The relevant length scale, which is the smaller of the screening radius (of the 
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photo-hole) and the escape depth, is ~ lOA. Thus t = 3 x lO"^'^ s, which should 
be compared with the time scale for electron-electron interactions (which are the 
dominant source of interactions at the high energies of interest): tee = {2tt)/u1p = 
4 X 10^^^ s, using a plasma frequency iOp — 1 eV for the cuprates (this would be 
even slower if c-axis plasmons are involved). If t <^ tee, then we can safely ignore 
vertex corrections. From our very crude estimate t/tee = 0.1, so that the situation 
with regard to the validity of the impulse approximation is not hopeless, but clearly, 
experimental checks are needed, and we present these in the next subsection. 

Very similar estimates can be made for renormalizations of the outgoing photo- 
electron due to its interaction with the medium; again electron-electron interactions 
dominate at the energies of interest. The relevant length scale here is the escape 
depth, which leads to a process of self-selection: those electrons that actually make 
it to the detector with an appreciable kinetic energy have suffered no collisions in 
the medium. Such estimates indicate that the "inelastic background" must be small 
and we will show how to experimentally obtain its precise dependence on k and uj 
later. 

Finally, one expects that final state linewidth corrections are small for quasi- 
2D materials based on the estimates made by Smith et al. |l^ In fact, this is yet 
another reason for the ease of interpreting ARPES data in quasi-2D materials. A 
clear experimental proof that these effects are negligible for Bi2212 will be presented 
later, where it will be seen that deep in the superconducting state, a resolution 
limited leading edge is obtained for the quasiparticle peak. 

2.4 Spectral Functions and Sum Rules 

Based on the arguments presented in the preceding subsection, we assume the valid- 
ity of the sudden approximation and ignore both the final-state linewidth broaden- 
ing and the additive extrinsic background. Then (B) is the only term that survives 
from all the terms described in Fig. ^ and the ARPES intensity is given by [ po| , pl[ 

/(k,^) =/o(k)/(c^)^(k,a;) (6) 

where k is the initial state momentum in the 2D Brillouin zone and uj the energy 
relative to the chemical potential. The prefactor /o(k) includes all the kinematical 
factors and the square of the dipole matrix element (shown in Eq. ^), /(tj) is 
the Fermi function, and A{]f.,uj) is the one-particle spectral function which will be 
described in detail below. 

We first describe some general consequences of Eq. ^ based on sum rules and 
their experimental checks. The success of this strategy employed by Randeria et 
al. j2^ greatly strengthens the case for a simple A(k,a;) interpretation of ARPES 
data. 

The one-particle spectral function represents the probability of adding or remov- 
ing a particle from the interacting many-body system, and is defined as j4(k,tj) = 
— (l/7r)lmG(k, a; -|- iO"'') in terms of the Green's function. It can be written as the 
sum of two pieces j4(k, ua) — (k, cj) + A+(k, ui), where the spectral weight to add 
an electron to the system is given by A+{]<.,uj) = X/m n e~'^^™K'^|Ck|fTi)p5(ct;-|- 
£m.—£n), and that to extract an electron is A_(k, tj) — Z^^ ^ e~''^'"|(n|ck|m)|^5(a;-|- 
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£n — £m)- Here |m) is an exact eigenstate of the many-body system with en- 
ergy £rn, Z is the partition function and /3 = 1/T. It follows from these defi- 
nitions that A_(k,aj) — f{u))A(k,u!) and A-|-(k,aj) = [1 ~ f{uj)] A(k,Lj), where 
/(oj) — l/[exp(/3aj) + 1] is the Fermi function. Since an ARPES experiment in- 
volves removing an electron from the system, the simple golden rule Eq. ^ can be 
generalized to yield an intensity proportional to A-(k,ij) 

We now discuss various sum rules for A(k,ij) and their possible relevance to 
ARPES intensity /(k,tj). While the prefactor Jo depends on k and also on the 
incident photon energy and polarization, it does not have any significant or T 
dependence. Thus the energy dependence of the ARPES lineshape and its T depen- 
dence are completely characterized by the spectral function and the Fermi factor. 
The simplest sum rule du)A(]i,uj) = 1 is not useful for ARPES since it involves 
both occupied (through A-) and unoccupied states {A+). Next, the density of states 
(DOS) sum rule '^^A(k,uj) — N{lo) is also not directly useful since the prefactor 
lo has very strong k-dependence. However it may be useful to k-sum ARPES data 
Np{u}) = 7o(k)^(k, oj) in an attempt to simulate the angle-integrated photoe- 
mission intensity. 

The important sum rule for ARPES is 



which directly relates the energy-integrated ARPES intensity to the momentum 
distribution n(k) — (cj^Ck). (The sum over spins is omitted for simplicity). Some- 
what surprisingly, the usefulness of this sum rule has been overlooked in the ARPES 
literature prior to Ref. pl[ ]. 

We first focus on the Fermi surface k — k^. One of the major issues, that we 
will return to several times in the remainder of this article, will be the question of 
how to define "k^" at finite temperatures in a strongly correlated system which 
may not even have well-defined quasiparticle excitations, and how to determine 
it experimentally. For now, we simply define the Fermi surface to be the locus of 
gapless excitations in k-space in the normal state, so that ^(kF,(x)) has a peak at 
w = 0. 

To make further progress with Eq. (Q), we need to make a weak particle- 
hole symmetry assumption: A(kF,—(^) ~ A(kF,i^) for "smaU" u, where "small' 
means those frequencies for which there is significant T-dependence in the spectral 
function. It then follows that |^] dn{kF)/dT — 0, i.e., the integrated area un- 
der the EDC at k^ is independent of temperature. To see this, rewrite Eq. ^ 
as n(kF) = \ - \ duJta.n\l{uJ/2T)[A{\^F,0J) - A{\iF,-uj)], and take its T- 
derivative. It should be emphasized that we cannot say anything about the value of 
n(kF), only that it is T-independent. (A much stronger assumption, A(}^F, —uj) = 
ACkFjU}) for all uj, is sufficient to give n(kF) = 1/2 independent of T). We empha- 
size the approximate nature of the kF-sum-ruIe since there is no exact symmetry 
that enforces it. 

We note that we did not make use of any properties of the spectral function 
other than the weak particle-hole symmetry assumption, and to the extent that 
this is also valid in the superconducting state, our conclusion dnCkp) / dT = holds 




(7) 
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equally well below Tc. There is the subtle issue of the meaning of "ki?" in the 
superconducting state. In analogy with the Fermi surface as the "locus of gapless 
excitations" above Tc, we can define the "minimum gap locus" below Tc. We will 
describe this in great detail in Sect. 5.1 below; it suffices to note here that "k^" 
is independent of temperature, within experimental errors, in both the normal and 
superconducting states of the systems studied thus far . 

In Fig. |(a) we show ARPES spectra for near-optimal Bi2212 {Tc = 87 K) at 
\^F along (vr, 0) to (tTjTt) at two temperatures: T = 13 K, which is well below Tc, 
and T = 95 K, which is in the normal state. The two data sets were normafized in 
the positive energy region, which after normalization was chosen to be the common 
zero baseline. (The essentially CLi-independent emission at positive energies, which is 
not cut-off by the Fermi function, is due to higher harmonics of the incident photon 
beam, called "second order light"). For details, see Ref. pl[ ]. 

In Fig. |(b) we plot the integrated intensity as a function of T and find that, in 
spite of the remarkable changes in the lineshape from 95K to 13K, the integrated 
intensity at k^ is very weakly T-dependent, verifying the sum rule dnCkp) / dT — 0. 
The error bars come from the normalization due to the low count rate in the lj > 
background. 

Let us discuss several other potential complications in testing the T-independence 
of the integrated intensity at kf. Note that matrix elements have no effect on this 
result, since they are T-independent. The same is true for a T-independent additive 
extrinsic background. In an actual experiment the observed intensity will involve 
convolution of Eq. ^ with the energy resolution and a sum over the k-values within 
the momentum window. While energy resolution is irrelevant to an integrated quan- 
tity, sharp k-resolution is of the essence. 

The T-independence of the integrated intensity is insensitive to the choice of 
the integration cutoff at negative u, provided it is chosen beyond the dip feature. 
It has been observed experimentally that the normalized EDCs at k^ have identi- 
cal intensities for all uj < — lOOmeV. This is quite reasonable, since we expect the 
spectral functions to be the same for energies much larger than the scale associ- 
ated with superconductivity. The fact that one has to go to lOOmeV (much larger 
than Tc ~ lOmeV) in order to satisfy the sum rule suggests that electron-electron 
interactions are involved in superconductivity. 

2.5 Analysis of ARPES Spectra: EDCs and MDCs 

In the preceding subsections we have presented evidence in favor of a simple spectral 
function interpretation of ARPES data on the quasi 2D high Tc cuprates. In the 
process, we saw that the ARPES intensity in Fig. M has a nontrivial lineshape which 
has a significant temperature dependence. We now introduce some of the basic ideas 
which will be used throughout the rest of the article to analyze and understand the 
ARPES lineshape. 

The one-electron Green's function can be generally written as G~^{]i,ijj) — 
G(^^(k,tj) — E{]i,Lo) where G'o(k,t<;) = l/[w — £k] is the Green's function of the 
noninteracting system, Sk is the "bare" dispersion, and the (complex) self-energy 
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E(]i,ijj) — S'CkjUj) + iS" (k, Lj) encapsulates the effects of all the many-body inter- 
actions. Then using its definition in terms of ImG, we obtain the general result 

^(k,.) = i , (8) 

We emphasize that this expression is entirely general, and does not make any as- 
sumptions about the validity of perturbation theory or of Fermi liquid theory. 

New electron energy analyzers, which measure the photoemitted intensity as 
a function of energy and momentum simultaneously, allow the direct visualization 
of the spectral function, as shown in Fig. ^, and have also suggested new ways of 
plotting and analyzing ARPES data. 

(a) In the traditional energy distribution curves (EDCs), the measured intensity 
7(k,tij) is plotted as a function of lu (binding energy) for a fixed value of k; and 

(b) In the new ji^] momentum distribution curves (MDCs), J(k,Lj) is plotted at 
fixed a; as a function of k. 

Until a few years ago, the only data available were in the form of EDCs, and even 
today this is the most useful way to analyze data corresponding to gapped states 
(superconducting and pseudogap phases). These analyses will be discussed in great 
detail in subsequent sections. We note here some of the issues in analyzing EDCs 
and then contrast them with the MDCs. First, note that the EDC lineshape is non- 
Lorentzian as a function of u) for two reasons. The trivial reason is the asymmetry 
introduced by the Fermi function /(w) which chops off the positive uj part of the 
spectral function. (We will discuss later on ways of eliminating the effect of the 
Fermi function). The more significant reason is that the self energy has non-trivial 
u dependence and this makes even the full yl(k,aj) non-Lorentzian in u as seen 
from Eq. |^. Thus one is usually forced to model the self energy and make fits to 
the EDCs. At this point one is further hampered by the lack of detailed knowledge 
of the additive extrinsic background which itself has tj-dependence. (Although, as 
we shall see, the MDC analysis give a new way of determining this background). 

The MDCs obtained from the new analyzers have certain advantages in studying 
gapless excitations near the Fermi surface |2^j2^,^. In an MDC the intensity is 
plotted as a function of k varying normal to the Fermi surface in the vicinity of 
a fixed \i.F{6), where 9 is the angle parametrizing the Fermi surface. For k near 
kp we may linearize the bare dispersion ek — v%{k — kp), where v'p{0) is the bare 
Fermi velocity. We will not explicitly show the 6 dependences of kp, Vp, or other 
quantities considered below, in order to simplify the notation. 

Next we make certain simplifying assumptions about the remaining k-dependences 
in the intensity 7(k,cj). We assume that: (i) the self-energy E is essentially inde- 
pendent of k normal to the Fermi surface, but can have arbitrary dependence on 
6 along the Fermi surface; and (ii) the prefactor /o(k) does not have significant k 
dependence over the range of interest. It is then easy to see from Eqs. ^ and ^that 
7(k,tj) plotted as function k (with fixed lu and 6) has the following lineshape. The 
MDC is a Lorentzian; 

(a) centered at k — kp + [lo ~ S' (Lj)]/vp; with 

(b) width (HWHM) Wm = \E"{lj)\/v%. 

Thus the MDC has a very simple lineshape, and its peak position gives the 
renormalized dispersion, while its width is proportional to the imaginary self en- 
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ergy. The consistency of the assumptions made in reaching this conclusion may be 
tested by simply checking whether the MDC lineshape is fit by a Lorentzian or not. 
Experimentally, excellent Lorentzian fits are invariably obtained (except when one 
is very near the bottom of the "band" or in a gapped state|26[). 

Finally, note that the external background in the case of MDCs is also very 
simple. One can fit the MDC (at each lu) to a Lorentzian plus a constant (at worst 
Lorentzian plus linear in k) background. From this one obtains the value of the 
external background including its uj dependence. Now this (^-dependent background 
can be subtracted off from the EDC also, if one wishes to. Note that estimating 
this background was not possible from an analysis of the EDCs alone. 

3 The Valence Band 

The basic unit common to all cuprates is the copper-oxide plane, Cu02- Some 
compounds have a tetragonal cell, a = b, such as the Tl compounds, but most have 
an orthorhombic cell, with a and b differing by as much as 3% in YBCO. There 
are two notations used in the literature for the reciprocal cell. The one used here, 
appropriate for Bi2212 and Bi2201, has F — M along the Cu — O bond direction, 
with M = (7r,0), and F — X{Y) along the diagonal, with Y = {iT,n). The other 
notation, appropriate to YBCO, has F — X{Y) along the Cu — O bond direction and 
F—S along the diagonal. This difference occurs because the orthorhombic distortion 
in one compound is rotated 45° with respect to the other. The main effect of the 
orthorhombicity in Bi2212 and Bi2201 is the superlattice modulation along the b 
axis, with Qsl parallel to F — Y . Except when refering to this modulation, we 
will assume tetragonal symmetry in our discussions. For a complete review of the 
electronic structure of the cuprates, see Ref. p^ ]. 

The Cu ions are four fold coordinated to planar oxygens. Apical (out of plane 
oxygens) exist in some structures (LSCO), but not in others. Either way, the apical 
bond distance is considerably longer than the planar one, so in all cases, the cubic 
point group symmetry of the Cu ions is lowered, leading to the highest energy Cu 
state having d^2_y2 symmetry. As the atomic 3d and 2p states are nearly degen- 
erate, a characteristic which distinguishes cuprates from other 3d transition metal 
oxides, the net result is a strong bonding-antibonding splitting of the Cu 3d^2_y2 
and O 2p a states, with all other states lying in between. In the stochiometric (un- 
doped) material, Cu is in a configuration, leading to the upper (antibonding) 
state being half filled. According to band theory, the system should be a metal. 
But in the undoped case, integer occupation of atomic orbitals is possible, and 
correlations due to the strong on-site Coulomb repulsion on the Cu sites leads to 
an insulating state. That is, the antibonding band "Mott-Hubbardizes" and splits 
into two, one completely filled (lower Hubbard band), the other completely empty 
(upper Hubbard band) [ psf . 

On the other hand, for dopings characteristic of the superconducting state, a 
large Fermi surface is observed by ARPES (as discussed in the next section). Thus, 
to a first approximation, the basic electronic structure in this doping range can be 
understood from simple band theory considerations. The simplest approximation is 
to consider the single Cu 3d^2_y2 and two O 2p (x,y) orbitals. The resulting secular 
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equation is 



Ced 2tsm{kxa) —2tsm{kya)\ 

2tsm{k^a) Ep (9) 

— 2tsin(fcj,a) / 

where e^j is the atomic 3d orbital energy, tp is the atomic 2p orbital energy, and t is 
the Cu-0 hopping integral. Diagonalization of this equation leads to a non-bonding 
eigenvalue E„b ~ fp, and 



E± = £li±±i±^(l^L^)2+4f2(sin2(fc,a)+sin2(fc,a)) (10) 
where + refers to antibonding, - to bonding. This dispersion is shown in Fig. |^a) 
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Fig. 7. (a) Simple three band estimate of the electronic structure of the Cu-O plane 
states; (b) EDC showing the whole valence band at the (7r,0) point; c) Intensity 
map of the whole valence band obtained by taking the second derivative of spectra 
such as the one in (b). The orbitals in (b) are based on the three band model, where 
black and white lobes corresponding to positive and negative wavefunctions. 



In Fig. [7|(b), we show an ARPES spectrum obtained at the (7r,0) point of the 
Brillouin zone for Bi2212. Three distinct features can be observed: the bonding 
state at roughly -6eV, the antibonding state near the Fermi energy, and the rest of 
the states in between. This rest consists of the non-bonding state mentioned above, 
as well as the remainder of the Cu 3d and O 2p orbitals, plus states originating from 
the other (non Cu-O) planes. It is difficult to identify all of these "non-bonding" 
states, as their close proximity and broadness causes them to overlap in energy. 
Perhaps surprisingly, the overall picture of the electronic structure of the valence 
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band has the structure that one would predict by the simple chemical arguments 
given above, as shown in Fig. |^c). The most important conclusion that one can 
derive from Fig |^ is the early prediction by Anderson |^^, namely that there is 
a single state relevant to transport and superconducting properties. This state, 
the antibonding state in Fig ^, is well separated from the rest of the states, and 
therefore any reasonable theoretical description of the physical properties of these 
novel materials should arise from this single state. 

Despite these simple considerations, correlation effects do play a major role, 
even in the doped state. The observed antibonding band width is about a factor 
of 2-3 narrower than that predicted by band theory |^l[. In fact, the correlation 
effects are the ones we are most interested in, as they give rise to the many unusual 
properties of the cuprates, including the superconducting state. 

As mentioned above, in the magnetic insulating regime, the picture is quite 
different. The antibonding band splits into two, with the chemical potential lying 
in the gap Upon doping with electrons or holes, the chemical potential would 
move up or down. One would then expect to observe small Fermi surfaces, hole 
pockets centered at (7r/2, 7r/2) with a volume equal to the hole doping, x. These can 
be generated at the mean field level by folding Fig. ^a) back into the magnetic zone 
with Q=(7r,7r), and turning on an interaction, ?7e//, between the two folded bands. 
This can be contrasted with the large Fermi surface of the unfolded case, typically 
centered at (tt, tt) with a volume 1 + a;. It is still an open question whether there is a 
continuous evolution between these two limits, or whether there is a discontinuous 
change at a metal-insulator transition point. In any case, a proper description of the 
electronic structure must take the strong electron-electron correlations into account, 
even in the superconducting regime. 



4 Normal State Dispersion and the Fermi Surface 

The Fermi surface is one of the central concepts in the theory of metals, with elec- 
tronic excitations near the Fermi surface dominating all the low energy properties of 
the system. In this Section we describe the use of ARPES to elucidate the electronic 
structure and the Fermi surface of the high Tc superconductors. 

It is important to discuss these results in detail because ARPES is the only 
experimental probe which has yielded useful information about the electronic struc- 
ture and the Fermi surface of the planar Cu-O states which are important for high 
Tc superconductivity. Traditional tools for studying the Fermi surface such as the 
deHaas-van Alphen effect have not yielded useful information about the cuprates, 
because of the need for very high magnetic fields, and possibly because of the lack 
of well defined quasiparticles. Other Fermi surface probes like positron annihilation 
are hampered by the fact that the positrons appear to preferentially probe spatial 
regions other than the Cu-O planes. 

The first issue facing us is: what do we mean by a Fermi surface in a system 
at high temperatures where there are no well-defined quasiparticles? (Recall that 
quasiparticles, if they exist, manifest themselves as sharp peaks in the one-electron 
spectral function whose width is less than their energy, and lead to a jump discon- 
tinuity in the momentum distribution at T = 0.) Clearly, the traditional T = 
definition of a Fermi surface defined by the jump discontinuity in n(k) is not useful 
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for the cuprates. First, the systems of interest are superconducting at low temper- 
atures. But even samples which have low Tc's have normal state peak widths at 
Ef which are an order of magnitude broader than the temperature[ ^3|j2^ ] . If, as 
indicated both by ARPES and transport, sharp quasiparticle excitations do not 
exist above Tc, there is no possibility of observing a thermally-smeared, resolution- 
broadened, discontinuity in n(k). 

It is an experimental fact that in the cuprates ARPES sees broad peaks which 
disperse as a function of momentum and go through the chemical potential at a 
reeisonably well-defined momentum. We can thus adopt a practical definition of the 
"Fermi surface" in these materials as "the locus of gapless excitations" . 

Historically, the first attempts to determine the Fermi surface in cuprates were 
made on YBCOjsij, however, surface effects as well as the presence of chains appear 
to complicate the picture, so we will focus principally on Bi2212 and Bi2201, which 
have been studied the most intensively. Other cuprates which have also been studied 
by ARPES, include the electron-doped material NCCO ^] and, more recently, 
LSCO as a function of hole doping |3q ]. 

We discuss below various methods used for the determination of the spectral 
function peaks in the vicinity of Ef. In addition, we supplement these methods with 
momentum distribution studies, taking due care of matrix element complications. 
We will then discuss three topics: the extended saddle-point in the dispersion, the 
search for bilayer splitting in Bi2212, and (in Section 6.4) the doping dependence 
of the Fermi surface. 



4.1 Normal State Dispersion in Bi2212: A First Look 

We begin with the results obtained by using the traditional method of deducing 
the dispersion and Fermi surface by studying the EDC peaks as a function of mo- 
mentum. This method was used for the cuprates by Campuzano et al.^^, Olson 
et al.\il\, and Shen and Dessau |^], culminating in the very detailed study of Ding 
et a/. 1 37 1 . The use of EDC peak dispersion has some limitations which we discuss 
below. Nevertheless, it has led to very considerable understanding of the overall elec- 
tronic structure, Fermi surface, and of superlattice effects in Bi2212, and therefore 
it is worthwhile to review its results first, before turning to more refined methods. 

The main results of Ding et al.^^ on the electronic dispersion and the Fermi 
surface in the normal state (T = 95K) of near-optimal OD Bi2212 (Tc = 87K) 
using incident photon energies of 19 and 22 eV are summarized in Fig. ^ The 
peak positions of the EDO's as a function of k are marked with various symbols in 
Fig. B(b). The filled circles are for odd initial states (relative to the corresponding 
mirror plane), open circles for even initial states, and triangles for data taken in a 
mixed geometry (i.e. the photon polarization was at 45° to the mirror plane). The 
Fermi surface crossings corresponding to these dispersing states are estimated from 
the k-point at which the EDO peak positions go through the chemical potential 
when extrapolated from the occupied side. The kp estimates are plotted as open 
symbols in Fig. ^(a). 

We use the following square lattice notation for the 2D Brillouin zone of Bi2212: 
EM is along the OuO bond direction, with E = (0,0), M = (tt, 0), X = (tt, -tt) 
and Y — (tt, tt) in units of 1/a* , where a* = 3.83A is the separation between near- 
neighbor planar Ou ions. (The orthorhombic a axis is along X and b axis along 
Y). 
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In addition to tiie symbols obtained from data in Fig. |8| there are also several 
curves which clarify the significance of all of the observed features. The thick curve 
in Fig. |^(b) is e(k), a six-parameter tight- binding fit jss) to the dispersion data in the 
y-quadrant; this represents the main Cu02 "band". It cannot be overemphasized 
that, although this dispersion looks very much like band theory (except for an 
overall renormalization of the bandwidth by a factor of 2 to 3) , the actual normal 
state lineshape is highly anomalous. As discussed in Section 7 below, there are no 
well defined quasiparticles in the normal state. 

The thin curves in Fig. ^b) are e(k±Q), obtained by shifting the main band fit 
by ±Q respectively, where Q = (0.2l7r, 0.2l7r) is the superlattice (SL) vector known 
from structural studies . We also have a few data points lying on a dashed curve 
e(k + K^) with = {n,-K); this is the "shadow band" discussed below. 

The thick curve in Fig. BHa) is the Fermi surface contour obtained from the main 
band fit, while the Fermi surfaces corresponding to the SL bands are the thin lines 
and that for the shadow band is dashed. It is very important to note that the shifted 
dispersion curves and Fermi surfaces provide an excellent description of the data 
points that do not lie on main band. We note that the main Fermi surface is a large 
hole-like barrel centered about the (tt, tt) point whose enclosed area corresponds to 
approximately 0.17 holes per planar Cu. One of the key questions is why only one 
CuO main band is found in Bi2212 which is a bilayer material with two CuO planes 
per unit cell. We postpone discussion of this important issue to end of this Section. 

The "shadow bands" seen above, were first observed by Aebi et al. ^] in 
ARPES experiments done in a mode similar to the MDCs by measuring as a func- 
tion of k the intensity J^^ du)f{u})A{'k,u)) integrated over a small range Sco near 
a; = 0. The physical origin of these "shadow bands" is not certain at the present 
time. They were predicted early on to arise from short ranged antiferromagnetic 
correlations ^]. In this case the effect should become stronger with underdoping 
toward the AFM insulator, for which there is little experimental evidence |^^. An 
alternative explanation is that the shadow bands are of structural origin: Bi2212 
has a face-centered orthorhombic cell with two inequivalent Cu sites per plane, 
which by itself could generate a (tt, tt) foldback. Interestingly, it has been recently 
observed that the shadow band intensity is maximal at optimal doping p3[ . 

We now turn to the effect of the superlattice (SL) on the ARPES spectra. This 
is very important, since a lack of understanding of these effects has led to much 
confusion regarding such basic issues as the Fermi surface topology (see below), 
and the anisotropy of the SC gap (see Section 5). The data strongly suggest [M 
that these "SL bands" arise due to diffraction of the outgoing photoelectron off the 
structural superlattice distortion (which lives primarily) on the Bi-0 layer, thus 
leading to "ghost" images of the electronic structure at ek±q. 

We emphasize the use of polarization selection rules (discussed in Section 2.1) 
for ascertaining various important points. First, we use them to carefully check the 
absence of a Fermi crossing for the main band along (0,0) — (tt, 0), i.e. FM. Note 
that the Fermi crossing that we do see near (tt, 0) along FM in Fig. ^(a) is clearly 
associated with a superlattice umklapp band, as seen both from the dispersion data 
in Fig. |^(b) and its polarization analysis. This Fermi crossing is only seen in the 
FM _L (odd) geometry both in our data and in earlier work |^^. Emission from 
the main d^2_y2 band, which is even about FM, is dipole forbidden, and one only 
observes a weak superlattice signal crossing Ef. (We will return below to newer 
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data at different incident photon energies where the possibihty of a Fermi crossing 
along FM is raised again). 

Second, we use polarization selection rules to disentangle the main and SL 
bands in the X-quadrant where the main and umklapp Fermi surfaces are very 
close together; see Fig. ^(a). The point is that FX (together with the z-axis) and, 
similarly FY , are mirror planes, and an initial state arising from an orbital which 
has d^2_y2 symmetry about a planar Cu-site is odd under reflection in these mirror 
planes. With the detector placed in the mirror plane the final state is even, and one 
expects a dipole-allowed transition when the photon polarization A is perpendicular 
to (odd about) the mirror plane, but no emission when the polarization is parallel 
to (even about) the mirror plane. While this selection rule is obeyed along FY it 
is violated along FX. In fact this apparent violation of selection rules in the X 
quadrant was a puzz ling feature of all previous studies ^ of Bi2212. It was first 
pointed out in Ref. pq|, and then experimentally verified in Ref. [^, that this 
"forbidden" rX\\ emission originates from the SL umklapps. We will come back to 
the FX\\ emission in the superconducting state below. 

4.2 Improved Methods for Fermi Surface Determination 

We now discuss more recently developed methodologies for Fermi surface determi- 
nation. The need for these improvements arises in part from the practical difficulty 
of determining precisely where in k-space a state goes through Ef. This problem is 
particularly severe in the vicinity of the (tt, 0) point of the zone where one is beset 
by the following complications in both Bi2201 and 2212: 

(1) The normal state ARPES peaks are very broad. This has important implica- 
tions about the (non-Fermi-liquid) nature of the the normal state, as discussed later 
(section 7). 

(2) The electronic dispersion near (7r,0) is anomalously flat ("extended saddle- 
point"). 

(3) In Bi2212 the superlattice structure complicates the interpretation of the data 
as described above. Fortunately this complication is greatly reduced in Pb-doped 
Bi2201 and 2212. 

(4) The flnal complication comes from the strong variation of the k-dependent 
ARPES matrix elements with incident photon energy. This makes the use of changes 
in absolute intensity as a function of k to estimate Fermi surface crossings highly 
questionable. 

Note that the first two points are intrinsic problems intimately related to the 
physics of high Tc superconductivity, the third is a material-specific problem, while 
the last is specific to the technique of ARPES. Nevertheless, all of these issues 
must be dealt with adequately before ARPES data on Bi2212 can be used to yield 
definitive results on the Fermi surface. 

Eliminating the Fermi function: Recall that the peak of the EDC as a 
function of uj corresponds to that of f{uj)A{\f.,uj), which in general does not coincide 
with the peak of the spectral function ^(k, lo). If one has a broad spectral function, 
which at kp is centered about u) = Q, then the peak of the EDC will be at < 
(positive binding energy), produced by the Fermi function f{uj) chopping off the 
peak of A, in addition to resolution effects. 

Since the goal is to study the dispersive peaks in j4(k,tj), rather than in the 
EDC, one must effectively eliminate the Fermi function from the observed intensity. 
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We present two ways of acliieving tliis goal, and illustrate it with data on Bi2201 
where it permits us to study the broad and weakly dispersive spectral peaks (points 
(1) and (2) above) near (7r,0) without the additional complication of the superlat- 
tice. 

One approach is simply to divide the data by the Fermi function; more accu- 
rately one divides the measured intensity I(kp,u)) by f{u>) * R{uj), the convolution 
of the Fermi function with the energy resolution. Although this does not rigorously 
give the spectral function (because of the convolution), this is a good approxima- 
tion in situations where the energy resolution is very sharp. An excellent example 
of such an approach can be seen in the work of Sato et al.^^ on Pb-doped Bi2201 
(Bii.8oPbo.38Sr2.oiCuOti-^i) which is overdoped with a Tc < iK. In Fig. ^ we show 
the raw data in the vicinity of M at various temperatures in the left panel and 
the corresponding data "divided by the Fermi function" in the right panel. The 
"divided" data permits one to clearly observe the dispersion well above Ef partic- 
ularly at high temperatures, and thereby identify the Fermi crossings with a great 
degree of confidence. 

An alternative approach |^ to eliminating the Fermi function is to symmetrize 
the data. For each k define the symmetrized intensity by 7sym(kF, tj) = /(kf, t<j) + 
I{k.F,—u}) = lo Asymi^F , i^) ■ It is easy to show that 7ayin('^) will exhibit a local 
minimum, or dip, at oj = for an occupied k state, while it will show a local max- 
imum at = for an unoccupied k state. In practice, then, the Fermi crossing k^ 
is determined as follows: All EDCs along a cut are symmetrized and k^ is identi- 
fied as the boundary in k-space between points where Tgym has a local maximum 
versus a local minimum at u = 0. As shown in Ref. |9|, these arguments work 
even in the presence of finite resolution effects. We note that this method, and the 
one presented above, for eliminating the Fermi function require a very accurate 
determination of the chemical potential (zero of binding energy). 

In Fig. |l^ we show the results of a symmetrization analysis for an OD 
23K Bi2201 (Bii.6Pbo.4Sr2Cu06-a) sample. From the raw data along cuts parallel 
to (7r,0) to (7r,7r) (see Fig. 7 of Ref. g) and along (0,0) (27r, 0) (shown in 
Fig. 0(d)) one sees broad peaks whose dispersion is very flat near (7r,0), thus 
making it hard to determine kf from EDC dispersion alone. Nevertheless, the 
symmetrized data provide completely unambiguous results: in the top panels (a,b, 
and c) of Fig. |l^ we illustrate the use of symmetrized data to determine kf using the 
criterion described above. Two other features about this analysis are worth noting. 
First, on approaching kF from the occupied side, resolution effects are expected to 
lead to a flat topped symmetrized spectrum. Second, one expects an intensity drop 
in the symmetrized spectrum upon crossing k^, assuming that matrix elements are 
not strongly k-dependent. Both of these effects are indeed seen in the data and 
further help in deducing k^. 

It is equally important to be able to ascertain the absence of a Fermi crossing 
along a cut. In this respect, the raw data along (0, 0) — (tt, 0) — {2n, 0) in Fig. p^(d)is 
difficult to interpret: the "flat band" remains extremely close to Ef but does it cross 
Ef? It is simple to see from the symmetrized data in Fig. |lo|(e) the answer is "no". 
The symmetrized data do not show a peak centered at tj = for any k, thus 
establishing the absence of a Fermi crossing along this cut. 

Intensity Plots: A different method for determining the Fermi surface is to 
make a 2D plot as a function of k of J^^ dujl(k, lo), the observed intensity integrated 
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Fig. 9. 
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Fig. 10. (a,b,c) Symmetrized EDCs along cuts parallel to (tt, 0) to (tt, tt) for OD 
23K Bi2201 at T = 25K using hu = 22eV^ photons. The curve corresponding to 
kF along each cut is determined as explained in the text and is shown by a thick 
line. The lower panel (d) shows EDCs along the (0, 0) — > {2tt, 0) direction; panel (e) 
shows symmetrized curves corresponding to the data in (d) and shows the absence 
of a Fermi surface crossing along (0,0) to (7r,0). 



over a suitably chosen narrow energy interval Slu about the Fermi energy. At first 
sight this seems to be a very direct way to find out the k-space locus on which 
the low energy excitations live. However, as we discuss below, one has to be very 
careful in interpreting such plots since one is now focusing on the absolute intensity 
of the ARPES signal, which can be strongly affected by the k-dependence of the 
photoemission matrix elements. 

This approach was pioneered in the cuprates by Aebi and coworkers |^^, and 
in recent times with the availability of Scionta detectors with dense k- samplin g it 
has been used by several groups p^ , ^ , ^ , ^l|j5^ ] . We show as an example in Fig/O 
results from our group |5l| on an optimally doped Tc = 90K Bi2212 sample [psT 
using 33eV incident photons. 

Before commenting on the controversies about the Fermi surface crossing near 
the (tTjO) point, we first examine EDCs along the FY direction (middle panel of 
Fig. hil, where the left panel shows a two dimensional plot of the energy and momen- 
tum dependent intensity of photoelectrons along the FY cut). All of the features 
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- main band (MB), superlattice umklapp bands (UB) and the shadow band (SB) 

- seen before and discussed above are confirmed. We also see weaker, second 
order umklapps from the superlattice (corresponding to ±(0.42,0.42), twice the 
superlattice wavevector), which confirms the diffraction origin of the superlattice 
bands. 



i 

^ 



















Fig. 11. (a) Intensity /(k,(j;) and (b) EDCs along FY measured on an optimally 
doped sample (Tc=90K) at T=40K with 33 eV photons polarized along FX. (c) 
Integrated intensity (-100 to +100 meV) covering the X and Y quadrants of the 
Brillouin zone. Data were collected on a regular lattice of k points (spacing 1° along 
r'X and 0.26° along FY), (d) Integrated intensity (± 40 meV) as in (c), but in the 
normal state (T = 150K). Overlaid on (c) is the main band (black), ± umklapps 
(blue/red), and ± 2nd order umklapps (dashed blue/red lines) Fermi surfaces from 
a tight binding fit js^. 



We now turn to panels (c,d) of Fig. ^ where we plot the integrated inten- 
sity within a ±100 meV window about the chemical potential. We note the very 
rapid suppression of intensity beyond ~ O.S/'M jisj, which does not occur in data 
taken with 22eV incident photons. This has led some authors js^Js^ to suggest 
the existence of an electron- like Fermi surface with a crossing at this point. How- 
ever, Fretwell et al. and independently Borisenko et al. |^ have argued that this 
Fermi crossing along (0, 0) to (tt, 0) is actually due to one of the umklapp bands, 
and the near optimally doped Fermi surface is indeed hole like as earlier shown by 
Ding et aZ.||. 

This can be seen most clearly in Pb-doped samples, where the umklapp bands 
are not visible. In Fig. |l^, we show the Fermi energy intensity map of Borisenko et 
al. where the hole surface centered around (7r,7r) and its shadow band partner 
are quite apparent. 

One of the main reasons for the controversy surrounding the topology of the 
optimal doped Fermi surface is the fact that data taken at different incident photon 
energies hv lead to different intensity patterns. Our assertion, based on Refs. pl|,p[], 
is that the superlattice umklapp band is more noticeable at hi^ = 33eV compared 
with 22 eV since matrix element effects suppress the main band intensity at 33 
eV. For a detailed discussion about how to discriminate between a main band and 
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Fig. 12. Upper panel: Energy distribution map from the F— (vr, vr) direction in the 
BriUouin zone of Pb-doped BSCCO recorded at room temperature. Lower paneh 
momentum distribution map of Pb-doped BSCCO, recorded at room temperature 
(raw data). White horizontal dashed line represents a Icf-EDC, vertical ones cor- 
respond to an i?i^-MDC. In both cases the gray scale represents the photoemission 
intensity as indicated. The inset shows the three dimensional (kx,ky,uj) space, which 
is probed in ARPES of quasi-2D systems (from Ref . ) . 



superlattice Fermi crossing, we refer the reader to the cited papers, and also to 
Refs. 1^,^ for the use of polarization selection rules for this purpose. 

Matrix Element Effects: There is an important general lesson to be learned 
from the above discussion which is equally relevant for the n(k) methods to be 
discussed below. Changes in the ARPES intensities (either integrated over a small 
energy window or over a large energy range) as a function of k can be strongly 
affected by matrix element effects. This is, of course, obvious from the expression 
for the the ARPES intensity: /(k,tj) = /o(k; A)/(aj)^(k,a;). The key question is: 
after integration over the appropriate range in u), how do we differentiate between 
the k dependence coming from the matrix elements 7o (which we are not interested 
in per-se) from the k dependence coming from the spectral function? 

One possibility is to have a priori information about the matrix elements from 
electronic structure calculations But as we now show, even in the absence 

of such information, one can experimentally separate the effects of a strong k- 
variation of the matrix element from a true Fermi surface crossing. The basic idea 
is to exploit the fact that changing the incident photon energy one only changes the 
ARPES matrix elements and not the spectral function (or the resulting momentum 
distribution) of the initial states. 

We will use Bi2201 to illustrate our point since it has all the complications 
(points (1),(2) and (4)) listed above without the superlattice (point (3)). Fig. ^ 
(a), from Ref. shows the intensity (integrated over a large energy range) as 
a function of k for OD 23K Bi2201, and highlights the differences between data 
obtained at 22 eV and 34 eV incident photon energies. At 22 eV the maximum 
intensity occurs close to (vr, 0) and decreases both toward (0, 0) and (tt, tt), while at 
34 eV there is a strong depression of intensity on approaching (7r,0), resulting in a 
shift of the intensity maximum away from (tt, 0). 

From the discussion of Section 2.4, we can write the integrated intensity as: 
7(k) = /J*^ daj/(k, uj) — /o(k; u; A)7i(k). We attribute this loss in intensity around 
(7r,0) at 34 eV seen in Fig. |l3|(a) to strong k-dependence of Jo rather than n(k). 
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Experimentally we prove this by showing that the EDCs at the same point in the 
Brillouin zone obtained at the two different photon energies show exa ctly the same 
lineshape, i.e. one can be rescaled onto the other as shown in Fig. [L3(b). As an 
independent check of this, we have also shown that the symmetrization analysis 
leads to the same conclusion that there is no Fermi crossing along (0,0) — > (7r,0); 
see Fig. 9 of Ref. [|. 




Energj' (eV) 



Fig. 13. Bi2201-OD23K, (a) Integrated intensity (over the range is -500 to 100 
meV) along the (0, 0) (tt, 0) — > (tt, tt) directions for two incident photon energies 
hu =22 and 34 eV. (b) Comparison of the ARPES lineshape measured at 22 (dashed 
lines) and 34 eV (solid lines) at three different k points. One curve has been scaled 
by the multiplicative constant indicated to make it lie on top of the other. 

For completeness, we note that another possible source of incident /ii'-dependence 
in ARPES is fcz-dispersion. If there was c-axis dispersion, different photon energies 
would probe initial states with different fcz values consistent with energy conser- 
vation. However the scaling shown in Fig. |l^(b) proves that it is the same two- 
dimensional (fcz-independent) initial state which is being probed in the data shown 
here, and the /ij^-dependence arises entirely from the different final states that the 
matrix element couples to. 

Methods based on the Momentum Distribution Finally we turn to the 
use of the momentum distribution sum rule (discussed in Section 2.4) in deter- 
mining the Fermi surface. In principle, locating the rapid variation of w(k) offers a 
very direct probe of the Fermi surface which we emphasize is not restricted to Fermi 
liquids. (The T = momentum distribution for known non-Fermi liquid systems, 
such as Luttinger liquids in one dimension, do show a inflection point singularity 
at kf.) However in practice, one needs to be very careful about the k-dependence 
of the matrix elements, as clearly recognized in the original proposal . 

Here we discuss two approaches using the n(k) method to obtain informa- 
tion about the Fermi surface. The first method is to study the k-space gradient 
of the logarithm of the integrated intensity. The second method is to study the 
temperature-dependence of the integrated intensity and use the approximate sum 
rule ^n{\^F)/^T = discussed in Section 2.4. 



30 



J.C. Campuzano, M.R. Norman, and M. Randeria 



The gradient method was used in our early work [^,^ where Icf was estimated 
from the location of max |Vkn(k)|. The same method has also been successfully used 
later by other authors (MpTl. In the presence of strong matrix element effects, it is 
even more useful to plot the magnitude of the logarithmic gradient: | Vk/(k)|/7(k) 
which emphasizes the rapid changes in the integrated intensity. The logarithmic 
gradient filters out the less abrupt changes in the matrix elements and helps to 
focus on the intrinsic variations in nfkl . 

As an example we show in Fig. ^ the results |^ of such an analysis for an OD 
OK Bi2201 sample. In the top panels (a) and (c) we show the integrated intensity 
7(k) around the (7r,0) point obtained at two different photon energies: 22 eV and 
28 eV respectively. In the lower panels 11(b) and (d), we plot | Vk/(k)|//(k). Note 
that there are large differences between the two top panels, due to different matrix 
elements at 22 eV and 28 eV. However, as explained above, the logarithmic gradients 
in the bottom panels, which are more influenced by the intrinsic n(k), are much 
more similar. The Fermi surface can be clearly seen as two high intensity arcs 
curving away from the (tt, 0) point. Once matrix element effects are taken care of, 
the Fermi surface results obtained at the two different photon energies are quite 
similar, and in good agreement with the results obtained from independent methods 
like symmetrization on the same data set M. 
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Fig. 14. Bi2201-OD0K: (a and c) Integrated intensity (over -350 to +50 meV) J(k) 
measured at Tiv =22 and 28 eV around the (7r,0) point. Notice that the intensity 
maximum depends strongly upon the photon energy h,v.{h and d) Corresponding 
gradient of the logarithm, | Vk/(k)|/7(k), the maxima which correspond to Fermi 
crossings and clearly show that, independent of the photon energy, the Fermi surface 
consists of a hole barrel centered around (7r,7r). 
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Our final method for the determination of a Fermi crossing goes back to the sum 
rule that we had introduced earlier, dn{'kF)/dT — 0. Assuming that the matrix 
elements are T- independent on the temperature scales of interest, this immediately 
implies the integrated intensity at (and only at) is T- independent. In Section 
2.4, we had used this sum rule to get confidence in the validity of the single-particle 
spectral function interpretation of ARPES by verifying at assuming that ki? 
was known (by some other means). Now we can invert the logic: we can look at the 
T dependence of the integrated intensity, and identify kF as that point in k-space 
where the integrated intensity is T-independent. This is illustrated in Fig. ^ from 
the work of Sato et al. who determine k^ along the (tt, 0) — (tt, tt) cut for a 
highly OD Bi2201 sample. The result is completely consistent with that obtained 
by other methods, such as "division by the Fermi function" on the same sample 
(see Fig. |). 

EDCs vs MDCs: To conclude this discussion, let us note a very recent devel- 
opment for determining kp and the near-_E/ dispersion based on the MDCs, which 
are plots of the ARPES intensity as a function of k (in this case normal to the 
expected Fermi surface), at various fixed values of u. As shown in Section 2.5, the 
MDC peak position in the vicinity of the Fermi surface, i.e, near (k — lip,^ ~ 0) 
is given by: k = kp + [co ~ E'{lli)]/vp. 

Thus kp is determined by the peak location of the MDC at a; = 0. The fully 
renormalized Fermi velocity vf = Vp/[1 — dS' /duj] is given by the slope of the 
MDC peak dispersion. We note that the factor arising from the fc-dependence of 
the self-energy is already included in so that = ii^'"'''[l -|- dS' /de\i\. (To see 
this, note that the analysis of Section 2.5 can be easily generalized to retain the 
first order term dE' /de\i without spoiling the Lorentzian lineshape of the MDC 
provided this fc-dependence does not enter E"). 

As discussed earlier, the above results derive from the Lorentzian lineshape of 
the MDC which arises when three conditions are satisfied: the matrix elements do 
not depend on k, the self energy does not depend on k (except for the (fc — kp) 
variation of E' noted above) and dispersion can be linearized near the Fermi surface. 
The validity of these assumptions can be checked self-consistently by the Lorentzian 
MDC lineshape and the dispersion deduced from the data. 

The significance of this approach is that, as emphasized by Kaminski et aZ. j2^, 
the dispersions of the EDC and MDC peak positions are actually different in the 
cuprates; see Fig. ^(a) in Section 7.7. This difference arises due to the non-Fermi 
liquid nature of the normal state, so that the EDC peak dispersion is not given by 
the condition ui — v%{k ~ kp) ~ E' = Q but also involves in general E" . In contrast 
the MDC peak dispersion is rigorously described by the expression described above, 
and is much simpler to interpret. 

We expect that the MDC method for determining the dispersion and kp which 
has thus far been used mainly along the zone diagonal, will eventually be the 
method of choice, except when one is very close to the bottom of the band where 
linearization fails. 

4.3 Summary of Results on the Optimally Doped Fermi Surface: 

We have discussed a large number of methods for determination of the Fermi surface 
in Bi2201 and Bi2212 in the previous subsection. These include: (a) dispersion 
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of EDC peaks through Ef, (b) dispersion of peaks after division by the Fermi 
function, (c) symmetrization, (d) maps of intensities at Ef, (e) gradient of n(k), 
(f) T-dependence of n{k), and (g) MDC dispersion. In addition we also discussed 
using the /i!/-dependence of the data and polarization selection rules to eliminate 
matrix element effects and to identify superlattice Fermi crossings. 

The reader might well ask: why so many different methods? The reason is 
that the development of all of these methods has taken place to deal with the 
complications of accurately identifying the Fermi surface in the presence of the four 
problems listed at the beginning of the preceding subsection. Each method has its 
pros and cons, so that some, like (b) and (c) require very accurate Ef determination, 
which is not the case in (e) and (f) which use energy-integrated intensities. Most 
methods require dense sampling in k space, while method (f) requires in addition 
data at several temperatures. 

Given the complications of the problem at hand it is important to look for 
crosschecks and consistency between various ways of determining the Fermi surface. 
We believe that for optimally doped Bi2201 and 2212 there is unambiguous evidence 
for a single hole barrel centered about the (tt, tt) point enclosing a Luttinger volume 
of (1 + a;) holes where x is the hole doping. We discuss further below the issues of 
the doping-dependence of the Fermi surface and of bilayer splitting in Bi2212. 



4.4 Extended Saddle Point Singularity 

The very flat dispersion near the (tt, 0) point observed in all of the data is striking. 
Speciflcally, along (0, 0) to (tt, 0) there is an intense spectral peak corresponding 
to the main band, which disperses toward Ef but stays just below it. This is 
often called the "flat band" or "extended saddle point" , and appears to exist in all 
cuprates, though at different binding energies in different materials p^ , ^ ,^ . 

In our opinion this flat band is not a consequence of the bare electronic struc- 
ture, but rather a many-body effect, because a tight-binding description of such a 
dispersion requires fine-tuning (of the ratio of the next-near neighbour hopping to 
the near-neighbour hopping) which would be unnatural even in one material, let 
alone many. 

An important related issue is whether this flat band leads to a singular density 
of states. It is very important to recognize that, while Fig. |^(b) looks like a con- 
ventional band structure, the dispersing states whose "peak positions" are plotted 
are extremely broad, with a width comparable to binding energy, and these simply 
cannot be thought of as quasiparticles. This general point is true at all k's, but 
specifically for the fiat band region it has the effect of spreading out the spectral 
weight over such a broad energy range that any singularity in the DOS would be 
washed out. This is entirely consistent with the fact that other probes (tunneling, 
optics, etc.) do not find any evidence for a singular density of states either. 



4.5 Bilayer Splitting? 

On very general grounds, one expects that the two Cu02 layers in a unit cell of 
Bi2212 should hybridize to produce two electronic states which are even and odd 
under refiection in a mirror plane mid- way between the layers. Where are these two 



34 



J.C. Campuzano, M.R. Norman, and M. Randeria 



states? Why then did we find only one main "band" and only one Fermi surface in 
Bi2212? 

Let us first recall the predictions of electronic structure calculations |5^] . In sys- 
tems like Bi2212, the intra-bilayer hopping as a function of the in-plane momentum 
kis of the form p9|j60[ | f±(k) = — tz(cos fc^; — cos fc^)'^. Thus the two bilayer states are 
necessarily degenerate along the zone diagonal. However they should have a maxi- 
mum splitting at M = (tt, 0) of order 0.25 eV, which may be somewhat reduced by 
many-body interactions. 

Depending on the exact doping levels and on the presence of Bi-0 Fermi surface 
pockets, which are neither treated accurately in the theory nor observed in the 
ARPES data, we must obtain one of the two following situations: (1) the bilayer 
antibonding (A) state is unoccupied while the bonding (B) state is occupied at 
(tt, 0). This would lead to an A Fermi crossing along (0, 0) — (tt, 0) and a B Fermi 
crossing along (7r,0) — (7r,7r). As described at great length above we did not find 
evidence for a main band Fermi crossing along (0, 0) — (tt, 0) at least for the near 
optimal doped sample, therefore this possibility is ruled out. 

(2) The second possibility is that both the A and B bilayer states are occupied 
at the (7r,0). In this case, there should be two (in principle, distinct) Fermi cross- 
ings along (7r,0) — {n,n), although they might be difficult to resolve in practice. 
Nevertheless, one would definitely expect to see two distinct occupied states at the 
(tt, 0) point. Unfortunately, the normal state spectrum at (tt, 0) is so broad in the 
optimally doped and underdoped materials that it is hard to make a clear case for 
bilayer splitting. Thus an effort was made to search for this effect in the supercon- 
ducting state at T ^ Tc, when a sharp feature (quasiparticle peak) is seen (see 
Fig. hq) and one might hope that the bilayer splitting should be readily observable. 

The issue then is how to interpret the peak-dip-hump structure seen in the 
ARPES lineshape at (7r,0) in Fig. The peak-dip-hump structure will be dis- 
cussed at length in Section 7 below. Nevertheless, here we will briefiy address the 
question of whether: (I) the peak and the hump are the two bilayer split states, 
which are resolved below Tc once the peak becomes sharp? Or (II) is the non-trivial 
line shape due to many-body effects in a single spectral function A(k,ij)? 

Three pieces of evidence will be offered in favor of hypothesis (II) as opposed 
to (I), so that no bilayer spliting is observable even in the superconducting state of 
near optimal doped Bi2212. The first evidence comes from studying the polarization 
dependence of the ARPES matrix elements. For case (I) there are two independent 
matrix elements which, in general, should vary differently with photon polarization 
A, and thus the intensities of the two features should vary independently as A 
is varied. On the other hand, for case (II), the intensities of the two features are 
governed by a single matrix element. As shown in Fig. |l^ it was found in Ref. 
that by varying the z-component of A, the peak and hump intensities scale together, 
and thus the peak-dip-hump are all part of a single spectral function for Bi2212. 

A second piece of evidence comes from a comparison of the normal and 
superconducting state dispersions near the (vr, 0) point, which will be discussed in 
detail in connection with Figs. ^ and ^ of Section 7.4. From these data, we argue 
that there is no evidence for a feature above Tc, which would correspond to the 
dispersionless quasiparticle peak below Tc- Thus the dispersionless peak must be of 
many-body origin. 
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The third and final piece of evidence comes from both ARPES and SIS tunnel- 
ing. In the ARPES data j62j shown in Fig. ^ of Section 7.6 one sees the striking 
fact that while the energy scales of both the quasiparticle peak and the (tt, 0) hump 
increase with underdoping, their ratio is essentialy doping independent. Since the 
location of the peak is the (maximum) superconducting gap - as discussed in detail 
in Section 5 - the (tt, 0) hump energy scales with the gap. SIS tunneling data also 
finds the same correlation on a very wide range of materials (includiiig some which 
have a single CuO plane) whose gap energies vary by a factor of 30 [^. This pro- 
vides very strong evidence that both the peak and hump are related to many-body 
effects and not manifestations of bilayer spliting. 

We note that Anderson had predicted that many-body efi'ects within a single 
layer could destroy both the quasiparticles and the coherent bilayer splitting in the 
normal state. But why the splitting should not be visible in the superconducting 
state, where sharp quasiparticles do exist, is not clear from a theoretical point of 
view. 

Recently, the above picture has been challenged by a number of authors ||,|6|j67| 
What has become clear is that bilayer splitting is indeed present for heavily over- 
doped Bi2212 samples, and has been seen now by several groups, including our own. 
In Fig. ^ we show (a) the bilayer split Fermi surfaces and (b) the bilayer split 
EDCs observed by the Stanford group for a heavily overdoped (rc=65K) Bi2212 
sample. Note that the bilayer splitting can even be seen in the umklapp bands. 
The resulting dispersion is reproduced in Fig. ^ where one sees that the momen- 
tum dependence of the splitting follows that expected from electronic structure 
considerations p£| , |60| ]. 

How this effect evolves as a function of doping, though, is still controversial. 
In particular, if it is present for optimal doping, it is difficult to resolve. Moreover, 
some authors who advocate bilayer splitting still argue that the peak/dip/hump 
structure in the superconducting state is largely a many-body effect as we 
adovcate here. This is supported by the fact that the same lineshape is seen in 
trilayer Bi2223 [^^, where one would expect three features if layer splitting were 
causing these effects. And similar lineshapes have been seen by tunneling in single 
layer systems, in support of a many-body interpretation. The reader is refered to 
Ref. |2[, where the issue of bilayer splitting is discussed in greater detail than here. 



5 Superconducting Energy Gap 

In this Section, we will first establish how the superconducting (SC) gap manifests 
itself in ARPES spectra, and then directly map out its variation with k along the 
Fermi surface. This is the only available technique for measuring the momentum 
dependence of the energy gap, and complements phase-sensitive tests of the order 
parameter symmetry [g^]. Thus ARPES has played an important role jToj, 
in establishing the d-wave order parameter in the high Tc superconductors | |69[ . 
At the end of the Section, we will discuss the doping dependence of the SC gap 
and its anisotropy, and the implications of this study for various low temperature 
experiments like thermal conductivity and penetration depth. 
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5.1 Particle-Hole Mixing 

To set the stage for the experimental results it may be useful to recall particle-hole 
(p-h) mixing in the BCS framework (even though, as we shall see in Section 7, there 
are aspects of the data which are dominated by many body effects beyond weak 
coupling BCS theory). The BCS spectral function is given by 

^(k, uj) = ulr/-K{{uj -^Ekf + r') + vlr/-K({Lo + E^f^r^) (ii) 

where the coherence factors are ~ 1 — iik ~ ^ ^k/^^k) and F is a phe- 
nomenological linewidth. The normal state energy ek is measured from Ef and the 
Bogoliubov quasiparticle energy is -Ek = ^/ + |zl(k)p, where /i(k) is the gap 
function. Note that only the second term in Eq. hll with the «k-coefficient, would 
be expected to make a significant contribution to the EDCs at low temperatures. 

In the normal state above Tc, the peak of A{k,uj) is at w = tk as can be seen 
by setting Zi = in Eq. |l^. We would thus expect to see in ARPES a spectral peak 
which disperses through zero binding energy as k goes through k^ (the Fermi sur- 
face). In the superconducting state, the spectrum changes from ek to Ek; see Fig. 
As k approaches the Fermi surface the spectral peak shifts towards lower binding 
energy, but no longer crosses Ef. Precisely at k^ the peak is at ui = |Zi(kF)|, which 
is the closest it gets to Ef. This is the manifestation of the gap in ARPES. Further, 
as k goes beyond k_F, in the region of states which were unoccupied above Tc, the 
spectral peak disperses back, receding away from Ef, although with a decreasing 
intensity (see Eq. This is the signature of p-h mixing. 

Experimental evidence for particle-hole mixing in the SC state was first given 
in Ref. |^^. In Fig. we show normal and SC state spectra for Bi2212 for k's 
along the cut shown in the inset. In the normal state data in panel (b) we see 
the electronic state dispersing through Ef-. the k's go from occupied (top of panel) 
to unoccupied states (bottom of panel). The normal state dispersion is plotted as 
black dots in Fig. ^ (b). The k^ obtained from this dispersion is in agreement 
with that estimated from the |Vk?i(k)| analysis of the normal state data shown in 
Fig. |l|(a). 

We see from Fig. bd (a) that the SC state spectral peaks do not disperse through 
the chemical potential, rather they first approach tj = and then recede away from 
it. The difference between the normal and SC state dispersions is clearly shown in 
Fig. |l|(b). _ 

There are three important conclusions to be drawn from Fig. E]l (b). First, 
the bending back of the SC state spectrum for k beyond k^ is direct evidence for 
p-h mixing in the SC state. Second, the energy of closest approach to tJ = is 
related to the SC gap that has opened up at the FS, and a quantitative estimate 
of this gap will be described below. Third, the location of closest approach to 
OJ = ("minimum gap") coincides, within experimental uncertainties, with the kf 
obtained from analysis of normal state data. 

In fact by taking cuts in k-space which which are perpendicular to the normal 
state Fermi surface one can map out the "minimum gap locus" in the SC state, or for 
that matter in any gapped state (e.g., the pseudogap regime to be discussed in the 
following Section) . We emphasize that particle-hole mixing leads to the appearance 
of the "minimum gap locus" and this locus in the gapped state gives information 
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Fig. 21. (a) Integrated intensity versus k from the normal state data of Fig. |2C| (b) 
shown by black dots gives information about the momentum distribution n(k) . The 
derivative of the integrated intensity is shown by the black curve (arbitrary scale). 
The Fermi surface crossing k = is identified by the minimum in the derivative, 
(b) Normal state dispersion (closed circles) and SC state dispersion (open circles) 
obtained from EDC's of Fig. M. Note the back-bending of the SC state dispersion 
for k beyond kf which is a clear indication of particle-hole mixing. The SC state 
EDC peak position at k^ is an estimate of the SC gap at that point on the Fermi 
surface. 



about the underlying Fermi surface. (By this we mean the Fermi surface on which 
the SC state gap appears below Tc). In fact, the observation of p-h mixing in the 
ARPES spectra is a clear way of asserting that the gap seen by ARPES is due 
to superconductivity rather than of some other origin, e.g., charge- or spin-density 
wave formation. 



5.2 Quantitative Gap Estimates 

The first photoemission measurements of the SC gap in the cuprates was by Imer et 
al. using angle-integrated photoemission, and by Olson and coworkers using 
angle-resolved photoemission. The first identification of a large gap anisotropy con- 
sistent with d-wave pairing was made by Shen and coworkers . Ding et al. 
subsequently made quantitative fits to the SC state spectral function to study the 
gap anisotropy in detail. 

We now discuss the quantitative extraction of the gap at low temperatures 
(T < Tc) following Ding et al.^. In Fig. ||, we show the T = 13K EDCs for an 
87K Tc sample for various points on the main band FS in the F-quadrant. Each 
spectrum shown corresponds to the minimum observable gap along a set of k points 
normal to the FS, obtained from a dense sampling of k-space |74|. We used 22 eV 
photons in a FY _L polarization, with a 17 meV (FWHM) energy resolution, and a 
k-window of radius 0.0457r/a*. 

The simplest gap estimate is obtained from the mid-point shift of the lead- 
ing edge of Bi2212 relative to a good metal (here Pt) in electrical contact with 
the sample. This has no obvious quantitative validity, since the Bi2212 EDC is a 
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spectral function while the polycrystalline Pt spectrum (dashed curve in Fig. ^2|) 
is a weighted density of states whose leading edge is an energy- resolution limited 
Fermi function. We see that the shifts (open circles in Fig. KS) indicate a highly 
anisotropic gap which vanishes in the nodal directions, confirming earlier results by 
Shen et al. |7C[| . These results are qualitatively similar to one obtained from the fits 
described below. 

Next we turn to modeling p^ , [n| the SC state data in terms of spectral func- 
tions. It is important to ask how can we model the non-trivial line shape (with the 
dip-hump structure at high uj) in the absence of a detailed theory, and, second, how 
do we deal with the extrinsic background? We argue as follows: in the large gap 
region near (7r,0), we see a linewidth collapse for frequencies smaller than ~ 3 A 
upon cooling well below Tc- Thus for estimating the SC gap at the low temperature, 
it is sufficient to look at small frequencies, and to focus on the coherent piece of the 
spectral function with a resolution-limited leading edge. (Note this argument fails 
at higher temperatures, e.g., just below Tc). This coherent piece is modeled by the 
BCS spectral function, Eq. hll. 

The effects of experimental resolution are taken into account via 

r r + oo 

/(k,cj)=7o/ dk' duj'R{Lj-ij')f{ij')A{k',uj') (12) 

J 5k J — oo 



where R{uj) , the energy resolution, is a normalized Gaussian and (5k is the k- window 
of the analyzer. In so far as the fitting procedure is concerned, all of the incoherent 
part of the spectral function is lumped together with the experimental background 
into one function which is added to the / above. Since the gap is determined by 
fitting the resolution-limited leading edge of the EDC, its value is insensitive to 
this drastic simplification. To check this, we have made an independent set of fits 
to the small gap data where we do not use any background fitting function, and 
only try to match the leading edges, not the full spectrum. The two gap estimates 
are consistent within a meV. Once the insensitivity of the gap to the assumed 
background is established, there are only two free parameters in the fit at each k: 
the overall intensity 7o and the gap |Z\|; the dispersion ek is known from the normal 
state study, the small linewidth F is dominated by the resolution. 

The other important question is the justification for using a coherent spectral 
function to model the rather broad EDC along and near the zone diagonal. As 
far as the early data being discussed here is concerned, such a description is self- 
consistent [ p2|]7l| , though perhaps not unique, with the entire width of the EDC 
accounted for by the large dispersion (of about 60 meV within our k-window) along 
the zone diagonal. More recent data taken along (0, 0) to (tt, tt) with a momentum 
resolution of 5k ~ 0.0l7r/a* fully justifies this assumption by resolving coherent 
nodal quasiparticles in the SC state; see Section 7.7. 

The gaps extracted from fits to the spectra of Fig. ^ are shown as filled sym- 
bols in Fig. For a detailed discussion of the the error bars (both on the the gap 
value and on the Fermi surface angle), and also of sample-to-sample variations in 
the gap estimates, we refer the reader to Ref. [Q. The angular variation of the gap 
obtained from the fits is in excellent agreement with | cos{kx) — cos(fcy)| form. The 
ARPES experiment cannot of course measure the phase of the order parameter. 
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but this result is in strong support of d^-2_yi pairing |69| . Such an order param- 
eter arises naturally in theories with strong correlations and/or antiferromagnetic 
spin fluctuations ItsI. Moreover, the functional form of the anisotropy we find is 
consistent with electrons in the Cooper pair residing on neighboring Cu sites. That 
is, ARPES gives information on the spatial range of the pair interaction which is 
difficult to obtain from other techniques. 

Finally we comment on the temperature dependence of the gap. Unfortunately 
with increasing temperature the linewidth grows and a simple BCS-type modeling 
(valid for T ^ Tc) of just the coherent part of the spectral function is not possible. 
While T-dependence of the leading edge can certainly be used as a rough guide, this 
estimate is affected by both the gap and the linewidth (the diagonal self-energy). 
We will discuss methods for modeling the SC state data in Section 7. 

For completeness, we add a remark clarifying the earlier observation of two 
nodes in the X-quadrant |2^], and the related non-zero gap along FX in the rX\\ 
geometry It was realized soon afterwards that these observations were re- 

lated to gaps on the superlattice bands |^^, and not on the main band. To prove 
this experimentally, the X-quadrant gap has been studied in the FX _L geometry 
fn] ] and found to be consistent with F-quadrant d^2_y2 result described above. 

5.3 Doping Dependence of the SC Gap 

There are two important issues to be addressed about the doping dependence: How 
does the maximum gap change with doping? and how does the gap anisotropy 
evolve with doping? 

ARPES results on underdoped samples show that the maximum gap increases 
fr^ , !?^ with decreasing hole concentration. For a recent compilation of gap versus 
doping results in Bi2212, see Fig. ^ of Section 7.6. Identical results have also 
been obtained by tunneling spectroscopy [Q. This was at first quite unexpected 
since Tc decreases as one underdopes from optimality. The fact that 2A/kBTc is 
not constant with doping and can become an order of magnitude larger than its 
BCS weak-coupling value is very clear evidence that the SC phase transition on the 
underdoping side is qualitatively different from the BCS transition, a point that we 
will return to in the next Section on pseudogaps. 

We next turn to the question of the SC gap anisotropy as a function of doping. 
In this case it is the behavior of the gap in the vicinity of the node which is the 
most important since all low temperature T ^ Tc thermodynamic and transport 
properties are controlled by thermal excitation of quasiparticles in the vicinity of 
the nodes (where the SC gap vanishes on the Fermi surface). These low energy ex- 



citations have a Dirac-like spectrum J5(k) = \ v'^k\ -\- v'^k'f. where fcx (fcy) are the 



components of the k perpendicular (tangential) to the Fermi surface, and measured 
from the nodal point. Here vf is the nodal Fermi velocity controlling the dispersion 
perpendicular to the Fermi surface, while V2 = 9Zik/9fc|| is the the gap slope at the 
nodal point. 

The density of states for low energy excitations is then given hj N{lu) — 
2uj / {■nvFV2) for uj <^ A (in units where Ti = 1). This leads to characteristic tempera- 
ture dependences for various low temperature properties with coefficients essentially 
determined hy vf and V2. Specifically for T <^Tc, the specific heat goes like C{T) — 
CiT /{vfV2)\ the thermal conductivity goes like k{T)/T = C2 {vf/v2 + V2/vf) and 
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the superfluid density goes hke Ps{T) = ps(0) — c^ivp / V2)T . Here ci and C2 are 
known constants while C3 contains (apriori unknown) muhiplicative Fermi liquid 
parameters po| , ^ . 

Clearly ARPES has a unique ability to independently measure both vp, using 
the dispersion of the nodal quasiparticle, and V2, using the slope of the SC gap at the 
node. One can then make detailed comparison with various bulk measurements, and 
thereby test the validity of the description of the low teniperature properties of the 
SC state in terms of weakly interacting quasiparticles |80[|8l| , ^ . This is particularly 
important given the lack of sharply defined quasiparticles in the normal state of 
the cuprates. 

With these motivations in mind, a detailed measurement of the shape of the 
superconducting gap in Bi2212 as a function of doping was carried out by Mesot et 
al.^^. Although these measurements were carried out at a time when the energy 
and momentum resolutions were lower than those currently available, they still give 
useful information because the gap over the full range of angles over the irreducible 
zone was measured. Using the simple BCS spectral function fits (described above) 
and taking into account the measured dispersion and the known energy and mo- 
mentum resolutions, Mesot et aZ.found the results shown in Fig. The simple 
d-wave gap A — Aq cos(20) (Fig. is modified by the addition of the first har- 
monic zik = /\max[Bcos(2(?!)) -|- (1 - B) cos(6<;/))] , with < B < 1. Note that the 
cos(6(/f)) term in the Fermi surface harmonics can be shown to be closely related to 
the tight binding function cos{2kx) — cos{2ky), which represents next nearest neigh- 
bours interaction, just as cos(2</)) is closely related to the near neighbor pairing 
function cos(kx) — cos(fcy). From Fig. ^ we find that while the overdoped data are 
consistent with B ~ 1, the parameter B decreases as a function of underdoping. 

Before discussing the doping dependence of the results in detail, let us first look 
at the comparisons between ARPES and other probes near optimality. Low tem- 
perature specific heat data on Bi2212 is not available, but the thermal conductivity 
has been measured by Taillefer and coworkers on an optimal Tc — 89K sample. 
Now one can use the k/T formula given above with the coefficient C2 — k%n/3hd, 
where n/dis the stacking density of Cu02 planes. It is important to note that C2 is 
not renormalized by either vertex corrections or Fermi liquid parameters ^] . Thus 
one infers vf jv^ = 19 from the thermal transport data |8^], which is in remarkable 
agreement with the ARPES estimate of vp /v2 — 20 coming from the measured 
values of vf ~ 2.5 x 10^ cm/s and V2 ~ 1.2 x 10® cm/s for the near optimal Tc = 87K 
sample of Fig. (We note that in Ref. we used the notation va for one-half 
the gap slope at the node, which is related to V2 defined above by V2 = 2vA/hkF-) 

The comparison of ARPES results with the slope of the superfiuid density 
Ps ~ obtained from penetration depth measurements is more complicated, as 
discussed in more detail in Ref. ISSl , for two reasons. Experimentally, there seem to 
be discrepancies between the dX~ /dT results of various groups, and theoretically 
the slope of ps is renormalized by doping-dependent Fermi liquid parameters [^,0 
which are not known apriori. These parameters characterize the residual interac- 
tions between the nodal quasiparticles in the superconducting state. Thus, e.g., even 
at optimality the slope d\~'^/dT obtained using the ARPES estimate of vf/v2 and 
ignoring Fermi liquid renormalization is almost three times as large as the exper- 
imentally measured value of dX~^/dT in near optimal Bi2212 [^^. This indicates 
the importance of Fermi liquid renormalizations in order to make a quantitative 
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comparison. For more details on the doping dependence of these renormalizations, 
see Ref. |80|||||. 

FinaUy, let us return to the doping dependence of the results of Fig. In 
contrast to the maximum gap (at the (tt, 0) — (tt, tt) Fermi surface crossing) increas- 
ing as a function of underdoping, noted earlier, we see the gap slope at the node, 
which determines «2, decreasing with underdoping. This is a result of the doping 
dependence of the B parameter introduced above. 

This surprising observation raises several questions. First, could the flattening 
at the node be, in fact, evidence for a "Fermi arc", a line of gapless excitations, in 
the underdoped materials, especially since such arcs are seen above Tc (See Sect. 
6.3 for further discussion of Fermi arcs in underdoped materials). Given the error 
bars on gap estimates in Fig. |2^, it is impossible to rule out arcs in all the samples. 
Nevertheless, it is clear that there are samples (especially OD87K, UD80K and 
UD75K) where there is clear evidence in favor of a point node rather than an arc 
at low temperatures. Furthermore, it is very important to note that a linear T 
dependence of ps(T) at low temperature, for all doping levels, in clean samples 
gives independent evidence for point nodes, at least in YBCO | p(| . 

Second, is the change in gap anisotropy intrinsic, or related to impurity scattering^]? 
We can eliminate the latter explanation on two grounds. The maximum gap in- 
creases as the doping is reduced, opposite to what would be expected from pair 
breaking due to impurities. Also, impurity scattering is expected to lead to a char- 
acteristic "tail" to the leading edge |fe8|, for which there is no evidence in the 
observed spectra. 

Thus the flattening near the nodes with underdoping does appear to be an 
intrinsic feature which may be related to the increased importance of longer range 
pairing interactions as one approaches the insulator. It would be of great interest 
to study the details of the doping dependence of the gap anisotropy with the new 
Scienta detectors which have greatly improved energy and momentum resolution. 



6 Pseudogap 

In this Section we describe one of the most fascinating developments in the study 
of high Tc superconductors: the appearance of a pseudogap above Tc which is seen 
most prominently on the underdoped side of the cuprate phase diagram. Briefly 
the "pseudogap" phenomenon is the loss of low energy spectral weight in a win- 
dow of temperatures Tc < T < T"; see Fig. |25[ The pseudogap regime has been 
probed by many techniques like NMR, optics, transport, tunneling, /iSR and spe- 
cific heat; for reviews and references, see Refs. | |89| ^ po|] . ARPES, with its unique 
momentum-resolved capab ilities, has played a central role in elucidating the pseu- 
dogap phenomenon ]9l| , |92| , |93| , [l2y94 l- 

We will discuss in this Section ARPES results on the anisotropy of the pseu- 
dogap, its T-dependence, its doping dependence, and its efi'ect on the normal state 
Fermi surface. We emphasize that for the most part we will focus on the "low en- 
ergy" or leading edge pseudogap, and only mention ARPES evidence for the "high 
energy pseudogap" toward the end. We will conclude the Section with a summary 
of the constraints put by the ARPES data on various theoretical descriptions of 
the pseudogap. 
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6.1 Pseudogap near (tt, 0) 

In the underdoped materials, Tc is suppressed by lowering the carrier (hole) con- 
centration as shown in Fig. ^ In the samples used by our group | p^ , p2|j94) un- 
derdoping was achieved by adjusting the oxygen partial pressure during annealing 
the float-zone grown crystals. These crystals also have structural coherence lengths 
of at least 1,250A as seen from x-ray diffraction, and optically flat surfaces upon 
cleaving, similar to the slightly overdoped Tc samples discussed above. We denote 
the underdoped (UD) samples by their onset Tc. the 83K sample has a transition 
width of 2K and the highly underdoped 15K and lOK have transition widths > 5K. 
Other groups have also studied samples where underdoping was achieved by cation 
substitution [|l|j92|. 

We now contrast the remarkable properties of the underdoped samples with the 
near-optimal Bi2212 samples which we have been mainly focusing on thus far. We 
will first focus on the behavior near the (tt, 0) point where the most dramatic effects 
occur, and come back to the very interesting k-dependence later. In Fig. ^ |7^] we 
show the T-evolution of the ARPES spectrum at the (tt, 0) — > (tt, tt) Fermi crossing 
for an UD 83K sample. At sufficiently high temperature, the leading edge of the 
UD spectrum at k_F and the reference Pt spectrum coincide, but below a crossover 
temperature T* ~ 180K the leading edge midpoint of the spectrum shifts below 
the chemical potential. In Fig. ^ one can clearly see a loss of low energy spectral 
weight at 120K and 90K. It must be emphasized that this gap-like feature is seen 
in the normal (i.e., non-superconducting) state for Tc — 83K < T <T* = 180K. 

The doping dependence of the temperature T*, below which a leading-edge 
pseudogap appears near (7r,0), is shown in Fig. |2^. Remarkably T* increases with 
underdoping, in sharp contrast with Tc, but very similar to the low temperature 
SC gap, a point we will return to at the end of the Section. The region of the phase 
diagram between Tc and T* is called the pseudogap region. 

It is important to emphasize that our understanding of the lightly UD samples 
(e.g., the UD 83K sample) is the best among the UD materials. In such samples all 
three regimes - the SC state below Tc, the pseudogap regime between Tc and T* 
and the gapless "normal" regime above T* - can be studied in detail. In contrast, 
in the heavily UD samples (e.g. the UD lOK and UD 15K samples), not only is 
the SC transition broad, one also has such low Tc's and such high r*'s that only 
the pseudogap regime is experimentally accessible. Nevertheless, the results on the 
heavily underdoped samples appear to be a natural continuation of the weakly 
underdoped materials and the results (the trends of gap and T*) on the low Tc 
samples are in qualitative agreement with those obtained from other probes (see 
Ref. HJoll). 

The T-dependence of the leading-edge midpoint shift appears to be completely 
smooth through the SC transition Tc. In other words, the normal state pseudogap 
evolves smoothly into the SC gap below Tc. Nevertheless, there is a characteristic 
change in the lineshape in passing through Tc associated with the appearance of of a 
sharp feature below Tc in Fig. |2^. This can be identified as the coherent quasiparticle 
peak for T <^Tc. The existence of a SC state quasiparticle peak is quite remarkable 
given that the normal state spectra of UD materials are even broader than at 
optimality, and in fact become progressively broader with underdoping. In fact, the 
low temperature SC state spectra near (7r,0) in the UD systems (see Fig. By of 
Section 7.6) are in many ways quite similar to those at optimal doping, with the 
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Fig. 26. ARPES spectra along the (tt, 0) ^ (n, n) direction for an 83 K underdoped 
sample at various temperatures (solid curves). The thin curves in each panel are 
reference spectra from poly crystalline Pt used to accurately determine the zero of 
binding energy at each temperature. 



one crucial difference that the spectral weight in the coherent quasiparticle peak 
diminishes rapidly with underdoping ]9^ , p7| . 

6.2 Anisotropy of the Pseudogap 

We have already indicated that the pseudogap above Tc near the (tt, 0) point of 
the zone evolves smoothly through Tc into the large SC gap below Tc, and thus the 
two also have the same magnitude. Since the SC gap has the d-wave anisotropy 
(discussed in detail in the preceding Section), it is natural to ask: what is the 
k-dependence of the pseudogap above Tc? 

The first ARPES studies | |9l] , ^ , p3| showed that the pseudogap is also highly 
anisotropic and has a k-dependence which is very similar to that of the SC gap below 
Tc- Later work [Q further clarified the situation by showing that the anisotropy has 
a very interesting temperature dependence. We now describe these developments 
in turn. 

In Fig. ^ josj we plot the leading edge shifts for three samples at 14K: the 
slightly overdoped 87K and UD 83K samples are in their SC states while the UD 
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lOK sample is in the pseudogap regime. The gap estimate for each sample was made 
on the "minimum gap locus" (explained earlier in the context of the SC gap; see 
further below). The large error bars on the UD lOK sample come from the difficulty 
of accurately locating the midpoint of a broad spectrum. Also there is a flattening 
of the gap near the node, a feature that we discussed earlier for the SC gap in UD 
samples. 

The remarkable conclusion is that the normal state pseudogap has a very similar 
k-dependence and magnitude as the SC gap below Tc. 




FS ansle 

Fig. 27. Momentum dependence of the gap estimated from the leading-edge shift in 
samples with T^'s of 87K (slightly overdoped), 83K (UD) and lOK (UD), measured 
at 14K. For the sake of comparison between samples we made vertical offsets so 
that the shift at 45° is zero; the offsets are — 3meV for the 83K and -|-2meV for the 
lOK sample. The inset shows the Brillouin zone with the large Fermi surface. 



6.3 Fermi Arcs 

The T-dependence and anisotropy of the pseudogap was investigated in more detail 
in Ref. motivated by the following question. Normal metallic systems are char- 
acterized by a Fermi surface, and optimally doped cuprates are no different despite 
the absence of sharp quasiparticles (see Section ^) . On the underdoped side of the 
phase diagram, however, how does the opening of pseudogap affect the locus of low 
lying excitations in k-space? 

In Fig. ^ we show ARPES spectra for an UD 83K sample at three k points on 
the Fermi surface for various temperatures. The superconducting gap, as estimated 
by the position of the sample leading edge midpoint at low T, is seen to decrease 
as one moves from point a near (tt, 0) to b to c, closer to the diagonal (0, 0) — > 
(7r,7r) direction, consistent with a d^2_y2 order parameter. At each k point the 
quasiparticle peak disappears above Tc as T increases, with the pseudogap persisting 
well above Tc, as noted earlier. 

The striking feature which is apparent from Fig. ^ is that the pseudogap at 
different k points closes at different temperatures, with larger gaps persisting to 
higher T's. At point a, near (7r,0), there is a pseudogap at all T's below 180K, at 
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Fig. 28. (a,b,c): Spectra taken at three k points in tlie Y quadrant of tlie zone 
(siiown in (d)) for an 83K underdoped Bi2212 sample at various temperatures 
(solid curves). The dotted curves are reference spectra from polycrystalline Pt (in 
electrical contact with the sample) used to determine the chemical potential (zero 
binding energy). Note the closing of the spectral gap at different T for different 
k's, which is also apparent in the plot (e) of the midpoint of the leading edge 
of the spectra as a function of T. Panels (f) show a schematic illustration of the 
temperature evolution of the Fermi surface in underdoped cuprates. The d-wave 
node below Tc (top panel) becomes a gapless arc above Tc (middle panel) which 
expands with increasing T to form the full Fermi surface at T* (bottom panel). 

which the Bi2212 leading edge matches that of Pt. As discussed above, this defines 
T* above which the the largest pseudogap has vanished within the resolution of 
our experiment, and a closed contour of gapless excitations - a Fermi surface - is 
obtained. The surprise is that if we move along this Fermi surface to point b the 
sample leading edge matches Pt at 120K, which is smaller than T*. Continuing 
to point c, about halfway to the diagonal direction, we find that the Bi2212 and 
Pt leading edges match at an even lower temperature of 95K. In addition, spectra 
measured on the same sample along the Fermi contour near the (0, 0) (n, tv) line 
shows no gap at any T (even below Tc) consistent with d^2_y2 anisotropy. 

One simple way to quantify the behavior of the gap is to plot the midpoint of the 
leading edge of the spectrum; see Fig. Ea(e) . We note that a leading edge midpoint 
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at a negative binding energy, particularly for k point c, indicates the formation 
of a peak in the spectral function at oj = at high T. Further, we will say that 
the pseudogap has closed at a k point when the midpoint equals zero energy, in 
accordance with the discussion above. A clearer way of determining this will be 
presented below when we discuss the symmetrization method, but the results will 
be the same. 

From Fig. ^ we find that the pseudogap closes at point a at a T above 180K, 
at point b at 120 K, and at point c just below 95 K. If we now view these data as 
a function of decreasing T, the picture of Fig. Ea(f) clearly emerges. With decreas- 
ing T, the pseudogap first opens up near (7r,0) and progressively gaps out larger 
portions of the Fermi contour. Thus one obtains gapless arcs which shrink as T 
is lowered, eventually leading to the four point nodes of the d-wave SC gap. The 
existence of such arcs is apparent from the first ARPES work on the pseudogap 
||9l| ], where it was noted that the Fermi contours in the pseudogap phase did not 
extend all the way to the zone boundary (see Fig. 

Whether the arcs shrink to a point precisely at Tc or below Tc is not clear from 
the existing data. As discussed in the preceding Section, we do believe that arcs do 
not survive deep into the SC state where there is point node at T ^ Tc in clean 
samples, as also evidenced by the linear T drop in the superfluid density at low T. 

We next turn to a powerful visualization aid that makes these results very 
transparent. This is the symmetrization method introduced in Ref. |^], which ef- 
fectively eliminates the Fermi function / from ARPES data and permits us to focus 
directly on the spectral function A. Given ARPES data described by|2^ ~ 
7o/(aj)yl(k, Lj) with the sum over a small momentum window about the Fermi 
momentum Rf, we can generate the symmetrized spectrum + Making 
the reasonable assumption of particle-hole (p-h) symmetry for a small range of uj 
and Ck, we have A(ek,uj) = A(—ek, — for jtjj, jej less than few tens of meV. It then 
follows, using the identity /(— oj) = 1 — f{uj), that I{io) + I{—uj) = y^,^ IoA(k, u) 
which is true even after convolution with a (symmetric) energy resolution function; 
for details see the appendix of Ref. j^. The symmetrized spectrum coincides with 
the raw data for lj < -2.2Teff, where 4ATeff is the 10%-90% width of the Pt 
leading edge, which includes the effects of both temperature and resolution. Non- 
trivial information is obtained for the range < 2.2Tsff, which is then the scale 
on which p-h symmetry has to be valid. We have extensively checked this method, 
and studied in detail the errors introduced by incorrect determination of the chem- 
ical potential or of (which lead to spurious narrow features in the symmetrized 
spectra), and the effect of the small (1° radius) k-window of the experiment (which 
was found to be small). 

In Fig. ^ we show symmetrized data for the UD 83K underdoped sample 
corresponding to the raw data of Fig. ^ To emphasize that the symmetry is put 
in by hand, we show the ui > curve as a dotted line. At k point a near (tt, 0) the 
sharp quasiparticle peak disappears above Tc but a strong pseudogap suppression, 
on the same scale as the superconducting gap, persists all the way up to 180K (T*). 
Moving to panels b and c in Fig. ^ we again see pseudogap depressions on the scale 
of the superconducting gaps at those points, however the pseudogap fills up at lower 
temperatures: 120K at b and 95K at c. In panel c, moreover, a spectral peak at zero 
energy emerges as T is raised. All of the conclusions drawn from the raw data in 
Fig. ^ are immediately obvious from the simple symmetrization analysis of Fig. 
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There are many important issues related to these results that will be taken up 
in Section 7.3 where we describe modeling the electron self-energy in the pseudogap 
state. We will discuss there the remarkable T-dependent lineshape changes and the 
T-dependence of the gap itself. Here we simply note that, without any detailed 
modeling, the data|94| clearly show qualitative differences in the T-dependence at 
different k-points. Near the (n, 0) point the gap goes away with increasing temper- 
ature with the spectral weight filling-in, but no perceptible change in the gap scale 
with T. On the other hand, at Icf points halfway to the node, one sees a suppression 
of the gap scale with increasing temperature. 
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We conclude this discussion with a brief mention of the imphcations of our 
results. We believe that the unusual T-dependence of the pseudogap anisotropy 
will be a very important input in reconciling the different crossovers seen in the 
pseudogap regime by different probes. The point here is that each experiment is 
measuring a k-sum weighted with a different set of k-dependent matrix elements 
or kinematical factors (e.g., Fermi velocity). For instance, quantities which involve 
the Fermi velocity, like dc resistivity above Tc and the penetration depth below 
Tc (superfluid density), should be sensitive to the region near the zone diagonal, 
and would thus be affected by the behavior we see at k point c. Other types of 
measurements (e.g. specific heat and tunneling) are more "zone- averaged" and will 
have significant contributions from k points a and b as well, thus they should see a 
more pronounced pseudogap effect. Interestingly, other data we have indicate that 
the region in the Brillouin zone where behavior like k point c is seen shrinks as 
the doping is reduced, and thus appears to be correlated with the loss of superfluid 
density Pq]. Further, we speculate that the disconnected Fermi arcs should have 
a profound influence on magnetotransport given the lack of a continuous Fermi 
contour in momentum space. 

6.4 Evolution of the Fermi Surface with Doping 

We now discuss the doping dependence of the normal state Fermi surface on the 
underdoped side of the phase diagram. The first issue to face up to is: can the 
Fermi arcs described above be a manifestation of a Fermi surface with small closed 
contours centered about (7r/2,7r/2)? Such hole-pockets enclosing x holes (per pla- 
nar Cu) are suggested by some theories of lightly doped Mott insulators M] as 
alternatives to the large Fermi surfaces containing (1 -I- x) holes which would be 
consistent with the Luttinger counting. 

The T-dependence of the arcs is by itself evidence against their being part of 
a pocket Fermi surface. Nevertheless, if there were such small hole pockets then 
one should observe two features in the ARPES data: a closure of the Fermi arc on 
the other side of (7r/2,7r/2), which would be clear evidence for a "shadow band"- 
like dispersion ((tt, 7r)-foldback of the main band) in the UD samples. In a variety 
of UD samples we have carefully searched for both these features and found no 
evidence for either [Q. However this is a tricky issue, given the very broad spectra 
and possible materials problems in the highly UD samples. Nevertheless, given the 
available evidence, the gapless arcs that we observe |M| are simply an intermediate 
state in the smooth evolution of d-wave nodes into a full Fermi surface. This smooth 
evolution was carefully checked on an UD 83K sample where a dense mapping was 
done in k space at T = 90K, revealing only a small Fermi arc just above Tc. 

The other issue related to the Fermi surface is: what is its doping dependence 
above T* where the pseudogap effects are absent. While one can easily compare the 
near optimal and lightly UD Fermi surfaces, the rapid rise of T* with underdoping 
does not permit us to address this question. However, one can study the "minimum 
gap locus" in any gapped state, in close analogy with the manner in which this was 
defined in the SC state; see Section 5. 

There is also a more fundamental reason to study the "minimum gap locus" in 
the pseudogap regime. One wants to know whether the pseudogap is "tied" to the 
Fermi surface, or if it has some other characteristic momentum Q (unrelated to k/). 



Photoemission in the High Superconductors 



59 





(7t. 71) 



(0,0) 



1 2 3 4 5 6 

EDC Number 



. * . . . 

1 ^ A ^ 



(K, 0) 



Fig. 31. (a) Midpoint shifts (dots) and widths (diamonds) for an UD 83K sample 
at a photon energy of 22 eV at 90K for a cut shown by the open dots in (b). (b) 
Fermi surface at 160K (sohd triangles) and minimum gap locus at 90K (solid dots) . 
Notice that the two surfaces coincide within error bars. The error bars represent 
uncertainties of Fermi crossings as well as possible sample misalignment. The solid 
curve is a rigid band estimate of the Fermi surface. 

In Fig. ^(a) we follow the dispersion of an UD 83K sample in the pseudogap 
regime. Moving perpendicular to the (expected) Fermi surface from occupied to 
unoccupied states, one finds that that the dispersion first approaches the chemical 
potential and then recedes away from it. This locates a k-point on the minimum 
gap locus. For a lightly UD sample we find in Fig. 0(b) that this locus in the 
pseudogap regime coincides, within experimental error bars, with the Fermi surface 
determined above T* where there is no pseudogap. The pseudogap is thus tied to 
the Fermi surface in the same way the SC gap is, and is in contrast with, say, charge 
or spin density waves, which are tied to other characteristic Q vectors. 

In the more heavily underdoped samples, it is not possible to compare the 
minimum gap locus in the pseudogap state with the Fermi surface above T* , or the 
minimum gap locus below Tc, since the latter two are not measurable with T* too 
high and Tc too low. Nevertheless, if one assumes, by continuity, that the minimum 
gap locus in the pseudogap state gives information about the Fermi surface that 
got gapped out, then even for an highly UD sample one finds a large underlying 
Fermi surface, satisfying the Luttinger count of (1 + x) holes per planar Cu |42| as 
shown in Fig. 

The same conclusion has been recently reached by the Dresden group . Fig. ^ 
reproduces their Fermi energy intensity maps as a function of doping, where a large 
Fermi surface (plus its shadow band image) is always visible. They argue, though, 
that the volume is not quite 1+x, and they attribute this difference to the presence 
of bilayer splitting. 

6.5 Low Energy vs High Energy Pseudogaps 

In all of the preceding discussion we have focussed on the "low energy" or leading- 
edge pseudogap. It is important to point out that the phrase pseudogap is (some- 
what confusingly) also used to describe a higher energy feature, which we call the 
"high energy pseudogap". 

The presence of a high energy pseudogap was evident in the first ARPES work 
on the pseudogap, reproduced in Fig. 53 M] . As the doping is reduced from optimal 
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Fig. 32. Fermi surfaces of the 87K, 83K, and 15K samples. All surfaces enclose a 
large area consistent with the Luttinger count (see text). The solid lines are tight 
binding estimates of the Fermi surface at 18%, 13%, and 6% doping assuming rigid 
band behavior. 
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Fig. 33. Basal plane projection of the normal-state (300 K) Fermi surface of 
Bi(Pb)-2212 from high-resolution ARPES. The Ef intensity (normahzed to the 
signal at w = 0.3 eV) is shown in color. The Tc of each sample is indicated. The 
raw data cover half of the colored area of each map and have been rotated by 180 
around the F point to give a better k-space overview. The line dividing raw and 
rotated data runs almost vertically for the UD76K map and from top left to bottom 
right in all other maps. The sketch shows the Fermi surface for the OD69K data set 
as yellow barrel-like shapes defined by joining the maxima of fits to the normalized 
Ef MDCs. (from Ref. @) 
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doped, a gap opens up in a region around the (vr, 0) points of the zone. The energy 
of thisgap is significantly higher than the leading edge gap emphasized in later 
work (Mbsl . The resulting dispersion of this high energy feature looks reminiscent 
of what is expected for a spin density wave gap. As the doping is further reduced, 
an energy gap then opens up along the (tt, tt) direction, and the material becomes 
truly insulating. 




Brillmiin Zone Symmetry Points 



Fig. 34. (a),(b) Map peak centroids vs k for Bi2Sr2Cai_2,Dy2,Cu20g+a thin films 
and deoxygenated Bi2Sr2CaCu20g+a bulk samples, respectively, with various hole 
doping levels, (a) Filled oval, 1% Dy near optimal doping with Tc = 85 K; gray 
diamond, 10% Dy underdoped with Tc = 65 K; gray rectangular, 17.5% Dy under- 
doped with Tc = 25 K; triangle, 50% Dy insulator, (b) Filled oval, 600 air annealed 
slightly overdoped with Tc = 85 K; gray diamond, 550 argon annealed underdoped 
with Tc = 67 K (from Ref.|||l). 



In Fig. [|6^ we show the temperature dependence of the (tt, 0) spectrum for 
an UD 89K sample. Note that there is no coherent quasiparticle peak until the 
system is cooled below Tc, with only a broad incoherent spectrum observed for all 
T > Tc- The leading edge pseudogap which develops below T* is difficult to see on 
the energy scale of this figure; the midpoint shift at 135 K is 3 meV. However, a 
higher energy feature can easily be identified by a change in slope of the spectra as 
function of binding energy; this is also very clear in the data of Fig. ^ of Section 
7.6. At the highest temperature T = 247K this feature is just due to the Fermi 
function cutoff, but in the pseudogap regime, this feature actually represents the 
onset of loss of spectral weight on a high energy scale, and hence may be called the 
"high energy pseudogap" . It can also be seen from Fig. ^ that the energy scale of 
this feature is very similar to that of the well-known (tt, 0)-hump of the peak-dip- 
hump structure seen in the SC state. This connection will be discussed in detail 
in Section 7.6 where we also argue that the high energy pseudogap and the hump 
have similar dispersions |b2[. 



6.6 Origin of the Pseudogap? 

We conclude with a summary of ARPES results on the pseudogap and a brief 
discussion of its theoretical understanding. 
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As described above, the low-energy (leading edge) pseudogap has the following 
characteristics. 

• The magnitude of the pseudogap near (7r,0), i.e., the scale of which there is 
suppression of low energy spectral weight above Tc, is the same as the maximum 
SC gap at low temperatures. Further both have the same doping dependence. 

• There is a crossover temperature scale T* above which the full Fermi surface of 
gapless excitations is recovered. The pseudogap near (7r,0) appears below T*. 

• The normal state pseudogap evolves smoothly through Tc into the SC gap as 
a function of decreasing temperature. 

• The pseudogap is strongly anisotropic with k-dependence which resembles 
that of the d-wave SC gap. The anisotropy of the pseudogap seems to be T- 
dependent leading to the formation of disconnected Fermi arcs below T*. 

• The pseudogap is "tied" to the Fermi surface, i.e., the minimum gap locus 
in the pseudogap regime coincides with the Fermi surface above T* and the 
minimum gap locus deep in the SC state, at least in those samples where all 
three loci can be measured. 

The simplest theoretical explanation of the pseudogap, qualitatively consistent 
with the ARPES observations, is that it arises to due pairing fluctuations above Tc 



|100,30|. The SC gap increases with underdoping while Tc decreases. Thus in the 



underdoped regime Tc is not controlled by the destruction of the pairing amp litude 



as in conventional BCS theory, but rather by fluctuations of the phase [101] of the 
order parameter leading to the Uemura scaling Tc ~ ps [^^. Even though SC order 
is destroyed at Tc, the local pairing amplitude survives above Tc giving rise to the 
pseudogap features. A natural mechanism for such a pseudogap coming from spin 



pairing in a doped Mott insulator exists within t he R VB framework |102|, with 
possibility of additional chiral current fluctuations ||l03|. 



More recently the pairing origin of the pseudogap has been challenged. Some 



experiments [104] have been argued to suggest a non-pairing explanation with a 



competition between the pseudogap and the SC gap. A specif ic re alization of this 



scenario is the staggered flux or d-density wave mechanism ] 105 ] in which T* is 
actually a phase transition below which both time-reversal and translational in- 
variance are broken. A more subtle phase transition with only broken time-reversal 



has also been proposed ]106] as the origin of the pseudogap. 

Although a qualitative understanding of some of the characteristics of the pseu- 
dogap within the non-pairing scenarios is not clear at this time, these theories make 
sharp predictions abo ut br oken symmetries below T* which can be tested. A very 



recent ARPES study [107| of circular dichroism finds evidence in favor of broken 
time reversal, thus casting some doubt on the pairing fluctuation ideas. The last 
word has clearly not been said on this subject, either theoretically or experimen- 
tally, and the origin of the pseudogap remains one of the most important open 
questions in the field of high Tc superconductors. 



7 Photoemission Lineshapes and the Electron 
Self-Energy 

Under certain conditions, which were discussed in Section 2, ARPES measures the 
occupied part of the single particle spectral function, A{k., ui)f{ui), with A = ImG/n 
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where G is the Greens function. The latter can be expressed as G~ = — — 
I]{]i.,u;) where tk is the single-particle energy (defined by the kinetic energy and 
single-particle potential energy terms of the Hamiltonian) and E is the Dyson self- 
energy (i.e., everything else). Often, this form is associated with a perturbative 
expansion used to estimate S, but of course the expression is itself tautological. 
The purpose of writing G in this form is that it isolates all many-body effects in the 
function S. An advatange of ARPES is that one has the possibility of extracting 
S directly from the data, allowing comparison to various microscopic predictions 
for E. 

One of the more trivial examples of this is when one fits ARPES data to de- 
termine the superconducting gap, A. For instance, the work described in Section 
5 used a broadened form of BCS theory to fit the leading edge of the spectra. 
This is equivalent to 17 = —iF + AI./{lj -f ek + iP), F = describing standard BCS 
theory. The advantage of this procedure is the actual gap function. A, is extracted 
from the data, rather than ill defined quantities, such as the often utilized leading 
edge shift (midpoint of the leading edge) which is not the same as A because of 
lifetime and resolution effects. When this is done, a Zi^ is obtained which has rather 
spectacular agreement with that expected for a d-wave order parameter. Although 
ARPES contains no phase information of the order parameter, the linear behav- 
ior of Ak along the Fermi surface near the gap zero (node) implies a sign change. 
Morevoer, ARPES has the additional advantage of determining the shape of Ak in 
the Brillouin zone, which gives important information on the spatial range of the 
pairing interaction |Q, as also discussed in Section 5. 

Even when fitting data at low temperatures including energy and momentum 
resolution, a non-zero F is always needed. The origin of this residual F is still 
debated. It is larger than what is expected based on impurity scattering, and cer- 
tainly larger than that implied by various conductivity probes (thermal, microwave, 
and infrared). Although the transport scattering rate is different from ImE (and 
in particular, only Umklapp processes contribute to electrical conductivity), the 
discrepancy is still large enough to be noticable, even when taking into account 
the fact that in the simple approximation being employed here, F represents some 
average of ImE over a frequency range of order A. 

Although it has been suggested that the residual 7^ is due to surface inhomogene- 
ity e ffects (in particular, a distribution of A due to local oxygen inhomogeneities 
108 1), a more likely possibility is that it is the same effect which is seen in normal 
metals like FiFe2. In the latter case, it was convincingly argued that this was the 
expected final state lifetime contamination effect when attempting to extract E 
from ARPES spectra [10!:]. Although the latter is expected to vanish in the pure 
2D limit, even small 3D effects can lead to a noticable effect, since final state life- 
times are large. For instance, in simple models, its contribution to F is of order 
{vl/vl)Ff , where vl is the c-axis velocity of the initial state, that of the final 



state, and Ff is the final state lifetime [|^. Since Ff is typically of order 1 eV, then 
a velocity ratio of only 0.01 is sufficient to cause a residual F of 10 meV. 

With this as an introduction, in this section, we desire to take a more serious 
look at the issue of extracting E from the data. The most commonly employed 
strategy is to come up with some model for E, and then see how well it fits the 
data, as illustrated by the simple example above. We will discuss this approach in 
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more detail later. We start, though, with discussing an alternate approach which 
we have recently advocated. 



7.1 Self-Energy Extraction 

Let us first assume we know A. Given that, we can easily obtain E. A Kramers- 
Kronig transform of A will give us the real part of G 



ReG{Lo) = P I diu'^^^ (13) 
u)' — to 



where P denotes the principal part of the integral. Knowing now both ImG and 
ReG, then E can be directly read off from the definition of G. 



2 



{ReGY + {ImG) 
(ReGrfilmGy 

To obtain ReG using Eq. |l3[ we need to know A for all energies. From ARPES, 
though, we only know the product of A and /. (While unoccupied states can be 
studied by inverse photoemission, its resolution at present is too poor to be useful 
for our purposes). This is not a limitation if an occupied k-state is being analyzed 
and one can either ignore the unoccupied weight or use a simple extrapolation 
for it (except that only ReE + e is determined) . On the other hand, one is usually 
interested in k vectors near the Fermi surface. Therefore a key assumption will have 
to be made. We can implement our procedure if we make the assumption of particle- 
hole symmetry, A{ek,u!) — A{—ek, —ui), within the small k-window centered at k^. 
Then, A is obtained by exploiting the identity A{e]^,uj)f{uj) + A{—e]^, —Lj)f{—Lj) — 
j4(ek,w), which holds even in the presence of the energy resolution convolution. 
Note, this can only be invoked at k.F, and was used previously to remove the 
Fermi function from ARPES data j94| , where it was denoted as the symmetrization 
procedure (note that the "symmetrized" data will correspond to the raw data for 
ui <~ — 2.2fcr). Although the particle- hole symmetry assumption is reasonable 
for small \uj\ where it can be tested in the normal state by seeing whether the 
"symmetrized" spectrum has a maximum at the Fermi energy {Ef), it will almost 
certainly fail for sufficiently large co > 0. Nevertheless, since we only expect to 
derive E for ui < 0, then the unoccupied spectral weight will affect the result only 
in two ways. The first is through the sum rule J dLoA{ijj) — 1 which must be used 
to eliminate the intensity prefactor of the ARPES photocurrent. From Eq. we 
see that violation of the sum rule will simply rescale ImE, but not ReE due to the 
oj — e factor. Our normalization, though, is equivalent to assuming Ukp— 0.5, and 
thus does not involve "symmetrized" data. The second influence comes from the 
Kramers-Kronig transformation in Eq. which is a bigger problem. Fortunately, 
the contribution from large u' > 0, for which our assumption is least valid, is 
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suppressed by 1/{lo' — lu). Further, for k^, ek=0 and thus ReE is not plagued by 
an unknown constant. 

When using real data, it is sometimes desirable to filter the noise out of the 
data, as well as to deconvolve the energy resolution, before employing the above 
procedure. These details can be found in Ref. |110[. Moreover, it is assumed that 
any "background" contribution (see Section 2) has been subtracted from the data 
as well. 

In Fig. ||a, we show T=14K symmetrized data for a Tc=87K Bi2212 overdoped 
sample at the (7r,0) point |110|. We note the important differences in this super- 
conducting state spectrum, compared with the normal state one (which can be fit 
by a simple Lorentzian), due to the opening of the superconducting gap, with the 
appearance of a sharp quasiparticle peak displaced from Ep by the superconduct- 
ing gap, followed by a spectral dip, then by a "hump" at higher binding energies. 
The resulting S is shown in Fig. Bqb and c. At high binding energies, one obtains 
a constant ImE with a very large value (~ 300 meV). Near the spectral dip, ImE 
has a small peak followed by a sharp drop. 

Despite this sharp drop below 70 meV, ImE remains quite large at low frequen- 
cies. Then, below 20 meV, there is a narrow spike in ImE. This is the imaginary 
part of the BCS self-energy, /{yj + iO"''), which kills the normal state pole at 
a;=0. The resulting 1/uj divergence of the real part ReS, which creates new poles 
at ±Zi=32meV, is easily seen in Fig. This is followed by a strong peak in 
ReS near the spectral dip energy, which follows from the Kramers-Kronig trans- 
formation of the sharp drop in ImE. The strong peak in ReE explains why the 
low energy peak in A is so narrow despite the large value of ImE. The halfwidth 
of the spectral peak is given by 7^ = zImE where = 1 — dReE/du) {z is the 
quasiparticle residue). In the vicinity of the spectral peak, z^^ is large (~9), giving 
a _r of ~14 meV. We note, though, that F is still quite sizeable, and thus the peak 
is not resolution limited, as discussed above. 

We can contrast this result with that obtained in the pseudogap phase. In 
Fig. p7|a, we show T—95K symmetrized data from a rc=85K underdoped Bi2212 
sample at the (vr, 0) — (vr, tt) Fermi crossing. One again sees (Fig. p7|b) a peak in ImE 
at ij=0, but it is broadened relative to that of the superconducting state, and the 
corresponding divergence of ReE (Fig. |3^) is smeared out. Such behavior would 
be consistent with replacing the BCS self-energy A'^/{uj + iO'^) by Zi^/(tj -I- i-To), 
and can be motivated by considering the presence of pair fluctuations above Tc, 
as will be discussed further below. Note from Fig. |3^ that although the equation 
oj — ReE{uj) = is still satisfied at ~ A, zImE is so large that the spectral 
peak is strongly broadened in contrast to the sharp peak seen below Tc. Actually, 
to a good approximation, the spectral function is essentially the inverse of ImE in 
the range <~ 2A. We can also contrast this case with data taken above T*, the 
temperature at which the pseudogap "disappears". In that case, the spectrum is 
featureless, and the peak in ImE is strongly broadened. As the doping increases, 
this peak in ImE disappears. Further doping causes a depression in ImE to develop 
around a; = 0, indicating a crossover to more Fermi liquid like behavior. 
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7.2 Temperature Dependence of S 

We now turn to the rather controversial issue of how the spectrum at (tt, 0) varies 
as a function of temperature, that is, how one interpolates between Figs. ^ and 
§. 

Fig. 38. 
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In Fig. we show data talcen for an optimal doped (Tc=90K) Bi2212 sample 
111|. The leading edge of the spectral peak is determined by the superconducting 
gap, whose energy stays fairly fixed in temperature, and persists above Tc (the 
pseudogap). On the trailing edge, one sees a spectral dip, whose energy also remains 
fixed in temperature, and becomes filled in above Tc due to broadening of the 
trailing edge of the peak. 
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In Fig. p9|a, ImU is plotted for various temperatures. At low temperatures 
and energies, it is again characterized by a peak centered at zero energy due to 
the superconducting energy gap, and a "normal" part which can be treated as 
a constant plus an uj'^ term. A maximum in ImS occurs near the energy of the 
spectral dip. Beyond this, ImS has a large, nearly frequency independent, value. 
As the temperature is raised, the zero energy peak broadens, the constant term 
increases, and the co^ term goes away. 

In Fig. ^9|b, the quantity uj ~ ReS is plotted. At low temperatures and energies, 
there is a l/ui term due to the energy gap, and a "normal" part which is linear 
in Lu. As expected, the zero crossing is near the location of the spectral peak. 
Beyond this, there is a minimum near the specral dip energy, then the data are 
approximately linear again, but with a smaller slope than near the zero crossing. 
As the temperature is raised, the gap (l/w) term broadens out and the low energy 
linear in ui term decreases, paralleling the behavior discussed above for ImS. 

From Figs. ^ and we see that rather than the spectral peak decreasing in 
weight with increasing temperature, it disappears by broadening strongly in energy. 
This can be seen directly by inspecting Fig. BS, in that as the temperature increases, 
ImS (Fig. |3E|a) in the vicinity of the peak increases in magnitude with T, and z^^ 
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(roughly the slope in Fig. |39p near the zero crossing) decreases with T. In fact, it is 
the strong T variation of ImE and , and the fact that they operate in concert, 
which is responsible for the rapid variation in the effective width of the spectral 
peak with T. 

The above analysis is important in that it shows how coherence is lost in the 
system. It is apparent from Figs. ^ and ^ that once a temperature is reached where 
the spectral peak is no longer discernable in the data, the difference in behavior 
of the self-energy between low energies and high energies is lost. That is, once the 
spectral dip is filled in, the low and high energy behaviors have merged, and the 
sharp peak and broad hump at low temperature is simply replaced by a single broad 
peak (with a leading edge gap due to the pseudogap). This is consistent with the 
spectral peak simply losing its integrity as the temperature is raised. The analysis 
does not support a picture of a well defined quasiparticle peak whose weight simply 
disappears upon heating, as has been suggested by other authors p&j- 



7.3 Modeling E 

This behavior can be further quantified by fitting the self-energy for binding energies 
smaller than the dip energy to that expected for a superfluid Fermi liquid, and 
exploring the temperature depe ndence of the resulting parameters. The reader can 



find this analysis in Ref. [111]. Rather, we will discuss here a simpler analysis 



we performed where we contrasted the temperature dependence of overdoped and 



underdoped samples [112]. In this case, we chose to look at data at the (tt, 0) — (tt, tt) 
Fermi crossing (antinode). The dip/hump structure is considerably weaker here 
than at the (tt, 0) point, allowing us to concentrate on more general aspects of the 
spectra. 

We begin with the overdoped sample, where there is no strong pseudogap effect. 
The simplest self-energy which can describe the low energy data at all T is 

r(k, c^) = -iA + A'/[{uj + iO+) + e(k)]. (15) 



Here A is a single-particle scattering rate taken, for simplicity, to be an tj-independent 
constant. It is effectively an aver age of the (actual oj-dependent) E" over the fre- 



quency range of the fit. (In Ref. [Ill , we generalized this to a constant plus an uj 
term, and included the resulting linear uj contribution to E'). The second term is 
the BCS self-energy discussed previously. 

In Fig. we show symmetrized data for an o verd oped Tc=82K sample at the 



antinode together with the fits obtained as follows [112[. Using Eq. |15| we calculate 
the spectral function 

^A(k,c^) = r"(k,a;)/ [{lo - ek - ^'(k,^;))^ + ^"(k,^^)^] , (16) 

convolve it with the experimental resolution, and fit to symmetrized data. The fit 
is restricted to a range of ±45 meV given the small gap in the overdoped case and 
the sharpness of the quasiparticle peaks below A. We find that Eq. |l^ describes 
the low energy data quite well. 
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The T- variation of the fit parameters A and A are shown in Fig. [4f|a. A{T) 
decreases with T, and although small at Tc, it only vanishes above Tc, indicating 
the possibility of a weak pseudogap. This effect is sample dependent, in that several 
overdoped samples we have looked at, the gap vanishes closer to Tc. We caution 
that the error bars shown in Fig. pf| a are based on the RMS error of the fits, but 
do not take into account experimental errors in /i and lip. 

is found to be relatively T- independent in the normal state. Below Tc, 
we see that A decreases very rapidly, and can be perfectly fit to the form a + 
bT^. This rapid drop in linewidth leading to sharp quasiparticle peaks at low T, 
which can be seen directly in the ARPES data, is consistent with microwave and 
thermal conductivity measurements, and implies that electron-electron interactions 
are responsible for Fi. Note the clear break in Fi at Tc, despite the fact A has not 
quite vanished. We have seen similar behavior to that described above for a variety 
of overdoped samples at several k points. 

We next turn to the more interesting underdoped case. We find that near (vr, 0) 
the self-energy (jl^) cannot give an adequate description of the data, in that it 
does not properly describe the pseudogap and its unusual "filling in" above Tc. 
Theoretically, we cannot have a divergence in S(kF,u! = 0) in a state without 
broken symmetry. A simple modification of the BCS self-energy rectifies both these 
problems: 

^(k, to) = -iFi + A^/[uj + e(k) + iFo]. (17) 



The new term Fo{T) should be viewed as th e inverse pair lifetime. The theoretical 
motivation for Eq. is given in Ref. [112|. We stress that this three parameter 
form is again a minimal representation of the pseudogap self-energy. Since it is not 
obviously a unique representation, it is very important to see what one learns from 
the fits. 

In Fig. ^o|b, we show symmetrized data at the antinode for a Tc=83K under- 
doped sample. Below Tc we see quasiparticle peaks. Above Tc these peaks disappear 
but there is still a large suppression of spectral weight around ijj=0. As T is raised 
further, the pseudogap fills in (rather than closing) leading to a fiat spectrum at 
a temperature of T* (200K). The self-energy, Eq. uM, gives a good fit to the data. 
These fits were done below Tc over a larger energy range (±75 meV) than in the 
overdoped case because of the larger gap. The range above Tc was increased to ±85 
meV so as to properly describe the pseudogap depression. 

In Fig. ^b, we show the T-dependence of the fit parameters. We find a number 
of surprises. First, A is independent of T within error bars. Similar behavior has 
been inferred from specific heat and tunneling data. This T-independence is in 
contrast to the behavior of the overdoped 82K sample with almost identical Tc at 
the same k point. In addition, for the underdoped sample, the gap evolves smoothly 
through Tc. 

The single-particle scattering rate Fi{T) for the underdoped sample is found 
to be qualitatively similar to the overdoped case. It is consistent with being T- 
independent above Tc, but with a value over twice as large as the overdoped case 
(allowing A to vary above Tc does not improve the RMS error of the fits). Second, 
we see the same rapid decrease in Fi below Tc as in the overdoped case. Note again 
the clear break at Tc. 
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The most interesting result is FoiT). We find 71) = below Tc and proportional 
to T — Tc above. This behavior is robust, and is seen in all the fits that we have 
tried. Moreover, a non-zero Fq is needed above Tc to obtain a proper fit to the 
data (its efi^ect cannot be reproduced by varying the other parameters). The fact 
that this T-dependence is exactly what one expects of an inverse pair lifetime is 
a non-trivial check on the validity of the physics underlying Eq. ^ Further, we 
observe from Fig. ^ that T* corresponds to where ^(21) ^ J^o{T). This condition 
can be understood from the small lj expansion of Eq. |l7| 

The next important question is whether the T-dependence at the antinode 
described above exists at other kF-points. To answer this, we have looked at T- 
dependent data for a number of underdoped samples at two different k vectors. 
All data at the antinode give results similar to those for the 83K sample. However 
at the second k-point, about halfway between the antinode and the node along 
the (0, 0) — (tt, tt) direction, we see quite different behavior. We demonstrate this 
in Fig. kda where symmetrized data for a 77K underdoped sample are shown. For 
the antinode, one clearly sees the gap fill in above Tc, with little evidence for any 
T-dependence of the position of the spectral feature defining the gap edge, just as 
for the 83K sample. In contrast, at the second k point, the gap is clearly closing, 
indicating a strong T-dependence of A. Similar behavior is seen in other underdoped 
samples with Tc between 75 and 85K. 

In Fig. [ii^j, we show the T-dependence of A obtained from fits (over a range of 
±66 meV) at the second k point for the 77K sample. A is found to be strongly T- 
dependent, being roughly constant below Tc, then dropping smoothly to zero above. 
The strong T-dependence of A makes it difficult to unambiguously determine To 
from the fits at this k-point. On theoretical grounds, we expect that, here too, 
there is a non-zero To, and the closing of the pseudogap is again determined by 
A{T) ~ To(T), however this condition is satified by the rapid drop in A{T), rather 
than the rise in To(T). For completeness, we also show A{T) for this sample at the 
antinode, which has a similar behavior to that of the 83K sample. 

We see that these results give further evidence for the unusual k-dependences 
first noted in Ref. joij. Strong pairing correlations are seen over a very wide T- 
range near (7r,0), but these effects are less pronounced and persist over a smaller 
T-range as one moves closer to the zone diagonal. This is clearly tied to the strong 
k-dependence of the effective interaction and the unusual (anomalously broad and 
non-dispersive) nature of electronic states near (7r,0). 



7.4 Peak/Dip/Hump - Experiment 

We now turn to a detailed discussion of the peak/dip/hump lineshape. As mentioned 
above, a very broad normal state spectrum near the (tt, 0) point of the zone evolves 
quite rapidly for T < Tc into a narrow quasiparticle peak, followed at higher binding 
energies by a dip then a hump, t he la tter corresponding to where the spectrum 



recovers to i ts normal state value 113 1. Similar effects are observed in tunneling 
spectra |114|]. 

In FigTBs, we show spectra for a Tc =87K Bi2212 sample along F — M — Z, i.e., 
(0, 0) — (tt, 0) — (27r, 0), in (a) the normal state (105 K) and (b) the superconducting 
state (13 K), from which we note two striking features |^l|. First, we see that the 
low energy peak in the superconducting state persists over a large range in k-space. 
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even when the normal state spectra have dispersed away from the Fermi energy. 
Second, when the hump in the superconducting state disperses, it essentially follows 
that of the normal state spectrum. This is accompanied by a transfer of weight to 
the hump from the low frequency peak, which is fairly fixed in energy. The same 
phenomena are also seen along M to Y (Fig. ^3|c). We will argue that the unusual 
dispersion seen in the superconducting state of Fig. ^ is closely tied to the lineshape 
change discussed earlier (Fig. p^ . 

The simplest explanation of the superconducting state spectra would be the 
presence of two bands (e.g., due to bilayer splitting), one responsible f or th e peak 
and the other for the hump. However, this explanation is untenable ill . First, 
if the sharp peak were associated with a second band, then this band should also 
appear above T^. But there is no evidence for it in the normal state data. Second, 
if the peak and hump were from two different bands, then their intensities must be 
governed by different matrix elements. However, we found that the intensities 
of both features scaled together as the photon polarization was varied from in to 
out of plane, as if they were governed by a common matrix element (Section 4.5). 
These arguments suggest that the unusual lineshape and dispersion represent a 
single electronic state governed by non-trivial many-body effects, as assumed in the 
previous discussion (Figs. |3^^9[). For more overdoped materials, though, bilayer 
splitting should be taken into account, as discussed in Section 4.5. 

Under this assumption, the data are consistent with a strong reduction of the 
imaginary part of the self-energy (ImS) at low energies in the superconducting 
state (Fig. pq ). If the scattering is electron-electron like in nature, then ImE at 
frequencies smaller than ~ 3A will be suppressed due to the opening of the su- 
perconducting gap |ll6| . On closer inspection, though, a more interesting story 
emerges. First, from Figs. || and ||, we see that the superconducting and normal 
state data match beyond 90 meV. From 90 meV, the dip is quickly reached at 70 
meV, then one rises to the sharp peak. Notice that since the width of the peak is 
around 20 meV, then the change in behavior of the spectra (from hump, to dip, 
to the trailing edge of the peak) is occuring on the scale of the energy resolution. 
That means that the intrinsic dip must be quite sharp. This implies that the large 
ImS at high energies must drop to a small value over a narrow energy interval to 
be consistent with the data, i.e., there is essentially a step in ImE. In fact, the 
data are not only consistent with a step in ImE, but the depth of the dip is such 
that it is best fit by a peak in ImE at the dip energy, followed by a rapid drop to a 
small value. This behavior can again be seen from the independent analysis shown 
in Figs. ^ and 

What are the consequences of this behavior in ImE? If ImE has a sharp drop 
at Lj, then by Kramers-Kronig transformation, ReE will have a sharp peak at uj 
(Fig. ^). This peak can very simply explain the unusual dispersion shown in Fig. 
as it will cause a low energy quasiparticle pole to appear even if the normal state 
binding energy is large. The most transparent way to appreciate this result is to note 
that a sharp step in ImE is equivalent to the problem of an electron interacting with 
a sharp (dispersionless) mode, since in that case, the mode makes no contribution to 
ImE for energies below the mode energy, and then makes a constant contribution 
for e nergies above. This problem has been treated by Engelsberg and Schrieffer 
117 1, and extended to the superconducting state by Scalapino and coworkers |ll8| . 
The difference in our case is that since the effect only occurs below Tc, it is a 
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consequence of the opening of the superconducting gap in the electronic energy 
spectrum, and thus of a collective origin, rather than a phonon. 

To facilitate comparison to this classic work, in Fig. W4 we plot the position of the 
low energy peak and higher binding energy hump as a function of the energy of the 
single broad peak in the normal state. This plot has a striking resemblance to that 
predicted for electrons interacting with a sharp mode in the superconducting state, 
and one clearly sees the low energy pole which we associate with the peak in ReS. 
On general grounds, the flat dispersion of the low energy peak seen in Fig. ^ is a 
combination of two effects: (1) the peak in ReS, which provides an additional mass 
renormalization of the superconducting state relative to the normal state, and thus 
pushes spectral weight towards the Fermi energy, and (2) the superconducting gap, 
which pushes spectral weight away. This also explains the strong drop in intensity 
of the low energy peak as the higher binding energy hump disperses. 

An important feature of the data is the dispersionless nature of the sharp peak. 
The mode picture discussed above would imply a dispersion of the peak from Ak 
to uj = u!o + Ak as the normal state binding energy increases (where uq is the mode 
energy). However, this dispersion turns out to be weak. From the data of Fig. 
we infer an wo = l.SAmax, uiq being essentially the energy separation of the peak 
and dip. Since Ak is known to be of the d^2_y2 form, then Ak should go to zero as 
we disperse towards the F point. Therefore, the predicted dispersion is only from 
A 

max to l.2,A^ax (32 to 42 meV). 
Since the dip/hump structure is most apparent at the (tt, 0) points, it is natural 
to assume that it ha s so mething to do with Q = (tt, tt) scattering, as discussed by 
Shen and Schrieffer |119|. But here, we find a new effect. If one compares the data 
of Figs. p3| b and p3|c, one sees that a low energy peak also exists along (tt, 0) — (tt, tt) 
for approximately the same momentum range as the one from (7r,0) — (0,0). That 
is, if there is a peak for momentum p, one also exists for momentum p + Q. This 
can be understood, since the self-energy equations for p and p + Q will be strongly 
coupled if Q scattering is dominant. 



7.5 Mode Model 

For now, though, we ignore the complication of momentum dependence. The low- 
est order contribution to electron-electron scattering is represented by the Feynman 
diagram shown in the inset of Fig. ^ In the superconducting state, each internal 
line will be gapped by A. This implies that the scattering will be suppressed for 
\uj\ < 3A. This explains the presence of a sharp quasiparticle peak at low tem- 
peratures. What is not so obvious is whether this in addition explai ns th e strong 
spectral dip. Explicit calculations show only a weak dip-like feature |{l2C[|. T o un- 
derstand this in detail, we equate the bubble plus interaction lines (Fig. [45| inset) 
to an "a^F" as in standard strong-coupling literature. In a marginal Fermi liquid 
(MFL) at T=0, a^F{f2) is simply a constant in fi. The effect of the gap is to force 
a^F to zero for SI < 2A. The question then arises where the gapped weight goes. 
It could be distributed to higher energies, but in light of the above discussion, we 
might expect it to appear as a collective mode inside of the 2A gap. For instance, if 
the bubble represents spin fluctuations, a sharp mode will appear if the condition 
1 - ?7xo(q, n) ^0 Ls satisfied for il < 2 A. These three cases (MFL, gapped MFL, 
gapped MFL plus mode) are illustrated in Fig. HB. 
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Fig. 46. (a) ImE for MFL (solid line), gapped MFL (dotted line), gapped MFL 
plus mode (dashed line) , and simple d-wave model (dashed-dotted line) . Parameters 
are a=l, cJc=200meV, zi=30meV (0 for MFL), no=2A, and /o=30meV. (b) Spec- 
tral functions (times a Fermi function with T=14K) convolved with a resolution 
gaussian of a— 7.5 meV for these four cases (e=-34meV). 



S is easy to obtain analytically if we ignore the complication of the supercon- 
ducting density of states from the k — q line of Fig. ^ and just replace this by a 
step function at A. The resulting I mE fo r the gapped MFL and gapped MFL plus 



mode models are shown in Fig. mm |121| in comparison to the normal state MFL. 



Note that structure in a F at il appears in E &t \uj\ — il + A due to the gap in the 
k — q line. Moreover, the MFL plus mode is simply the normal state MFL cut-off 
at 3 A (this is obtained under the assumption that all the gapped weight in a^F 
shows up in the mode). In contrast, the gapped MFL decays linearly to zero at 3 A. 
The Nambu spectral function is given by 

with (a complex) Z{uj) = 1 — E{lj)/uj. These are shown in Fig. ||b and were 
convolved with a gaussian of a— 7.5 meV, typical of high resolution ARPES, with 
a constant ImS (Fo) added for > Zi to reduce the size of the quasiparticle 
peak. We note that there is no dip as such for the gapped MFL model, whereas the 
addition of the mode causes a significant dip. The latter behavior is consistent with 
experiment. Moreover, the mode model has the additional advantage that ImS 
recovers back to the normal state value by 3A, which is also in agreement with 
experiment in that the normal and superconducting state spectra agree beyond 90 
meV (Figs. |6|and||). 

We contrast this behavior with that expected for a simple d-wave model. To 
a first approximation, this can be obtained by replacing the step drop in ImS in 



the MFL plus mode model with {\uj\ - A)^ for < 3Z\ [h2d]. This is shown in 
Fig. p^ a as well, with the resulting spectrum in Fig. p6|b. Only a weak dip appears. 
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Moreover, we have analyzed models with the exponent 3 replaced by some n and 
have found that n must be large to obtain a dip as strong as seen in experiment. 
Therefore, the upshot is that at the least, something similar to a step is required 
in IraE to be consistent with experiment. 

In principle, we could take the above MFL plus mode model and fit experiment 
with it. We consider a simpler model. There are several reasons for this. First, the 
MFL model has a number of adjustable parameters associated with it. There is the 
coupling constant (a), the cut-off frequency {toa), and the mode energy (which is 
not in general 2A). Moreover, the spectrum for k points near the (7r,0) point does 
not appear to be MFL-like in nature. We have found that the normal state Bi2212 
spectrum is fit very well by a Lorentzian plus a constant background in an energy 
range less than 0.5eV. This is also true for Bi2201 spectrum where the normal state 
can be accessed to much lower temperatures. 

In the resulting Lorentzian model, the normal state E is purely an imaginary 
constant, and F is a mode at zero energy. In the superconducting state, this mode 
gets pushed back to some energy within 2A. This model is artificial in the sense that 
all the self-energy is being generated by the mode. That is why we went through the 
above discussion motivating the mode more properly as a rearrangement of a^F 
due to the superconducting gap. In practice, though, the results are very similar 
to the MFL plus mode model, and has the further advantage of having the several 
parameters of that model collapse to just the mode strength (A ) and mode position 
(i7o) of the Lorentzian model. Moreover, analytic results can still be obtained for 
E when the superconducti ng de nsity of states for the k — q line of Fig. ^ is taken 
into account. The result is 121] 



ImE{i^) = rQN{\u\) + riN{\ui\ - flo), \u\ > + A 
= roN{\u;\), A<\lu\ < Ho + A 
= 0, \lu\ < A 



(19) 



where N{uj) = ui / \/lo^~~A^ is the BCS density of states, and 



7vReE{uj) = FoNi-Lj) In 
-f riiV(i7o - LJ) In 



^/Z^^A^l/A 



no-iJ+ ^y{iu - Oo)^ - A'2\/A 



(20) 



where it has again been assumed that A is a. real constant in frequency. An s-wave 
density of states has been used to obtain an analytic result. A d-wave density of 
states will not be that different. The advantage of an analytic result is that it is 
useful when having to take spectra and convolve with resolution to compare to 
experiment. Our results are not very sensitive to Fo, included again to damp the 
quasiparticle peak. 

The resulting real (Eq. pc| ) and imaginary (Eq. |l^ parts of E at (tt, 0) are 
shown in Fig. WM. Note the singular behaviors at A (peak energy) due to the Fo 
term and at Qo + A (dip energy) due to the A term. In both cases, step drops in 
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ImE would also give singularities in ReS. The advantage of peaks in ImE (due to 
the SC density of states) is that it makes the dip deeper in better agreement with 
experiment. In Fig. [irf j and ^c, we show a comparison of the resulting spectral 
function (convolved with the experimental energy and momentum resolution) to 
experimental data at (tt, 0) for both wide and narrow energy scans, where a step 
edge background with a gap of A is added to the calculated spectrum. The resulting 
agreement is excellent. 

It is interesting to note that the mode energy we infer from the data is 41 meV, 
equi valent to a magnetic resonant mode energy observed in YBCO [123 and Bi2212 
1 124 1 by neutron scattering data at Q = (7r,7r)- The Q dependence of this mode 
correlates well with the observations of Fig. ^ To explore this in greater detail, 
we now consider the doping dependence of the peak/dip/hump structure. 



7.6 Doping Dependence 

We show data along (tt, 0) (tt, tt) for an underdoped 75K sample in the super- 
conducting state (Fig. p^) and in the pseudogap state (Fig. [isfa) ]6^ . Below Tc, 
the sharp peak at low energy is essentially dispersionless, while the higher energy 
hump rapidly disperses from the (vr, 0) point towards the (tt, 0) — » (tt, tt) Fermi 
crossing seen above T* . Beyond this, the intensity drops dramatically, but there 
is clear evidence that the hump disperses back to higher energy. In the pseudogap 
state, the high energy feature also shows strong dispersion, much like the hump 
below Tc, even though the leading edge is non-dispersive like the sharp peak in the 
superconducting state. 

In Fig. ^ we show the dispersion of the sharp peak and hump (below Tc), 
for a variety of doping levels, in the vicinity of the (tt, 0) point along the two 
principal axes. The sharp peak at low energies is seen to be essentially non-dispersive 
along both directions for all doping levels, while the hump shows very interesting 
dispersion. Along (7r,0) (0,0) (Fig. p9|a), the hump exhibits a maximum, with 
an eventual dispersion away from the Fermi energy, becoming rapidly equivalent 
to the binding energy of the broad peak in the normal state as one moves away 
from the region near (tt, 0). In the orthogonal direction (Fig. p9|b), since the hump 
initially disperses towards the (tt, 0) — ^ (tt, tt) Fermi crossing, which is known to be 
a weak function of doping, one obtains the rather dramatic effect that the dispersion 
becomes stronger with underdoping. We also note that there is an energy separation 
between the peak and the hump due to the spectral dip. In essence, the hump 
disperses towards the spectral dip, but cannot cross it, with its weight dropping 
strongly as the dip energy is approached. Beyond this point, one sees evidence of 
the dispersion bending back to higher binding energy for more underdoped samples. 

Fig. [s^ a shows the evolution of the low temperature spectra at the (tt, 0) point 
as a function of doping. The sharp quasiparticle peak moves to higher energy, 
indicating that the gap increases with underdoping (although this is difficult to see 
on the scale of Fig. pda). We see that the hump moves rapidly to higher energy 
with underdoping. These trends can be seen very clearly in Fig. |5C|b, where the 
energy of the peak and hump are shown as a function of doping for a large number 
of samples. Finally, we observe that the quasiparticle peak loses spectral weight 
with increasing underdoping, as expected for a doped Mott insulator; in addition 
the hump also loses spectral weight though less rapidly. 
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Fig. 50. Dependence of energy scale on carrier density: (a) Doping dependence of 
the spectra (T=15K) at the (7r,0) point. The inset shows Tc vs. doping, (b) Doping 
dependence of T* , and the peak and hump binding energies in the superconducting 
state along with their ratio (c), as a function of doping, x. The empirical relation 
between and x is given by T^/T^"'' = 1 - 82.6(a; - 0.16)^ with r™''^=95K. For 
T* , solid squares represent lower bounds. 




^ Hole doping Igwl » 

Fig. 51. (a) Doping dependence of the low-T (14 K) coherent weight (za)- The 
dashed line is a guideline showing that za increases linearly on the underdoped side, 
and tapers off on the overdoped side, (from Ref. jo^) (b) The doping dependence 
of the superconducting peak ratio (SPR) is plotted over a typical Bi2212 phase dia- 
gram. The solid line is a guide to the eye. Horizontal error bars denote uncertainty 
in determining the doping level (±0.01); vertical error bars denote uncertainty in 
determining the SPR(±1.5%). AF, antiferromagnetic regime; SC, superconducting 
regime, (from Ref. ||9q ]) 
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This effect has recently been quantified in greater detail, where it was found 
that the spectral weight of the peak varies linearly with doping, as reproduced in 
Fig. 1^,^ . We remark that Ding et al. ^] also found the unusual relation that 
the product of the peak weight times the peak energy is constant with doping. 

The hump below Tc is clearly related to the superconducting gap, given the 
weak doping dependence of the ratio between the hump and quasiparticle peak 
positions at (vr, 0), shown in Fig. po|c. Tunneling data find this same correlation on 
a wide variety of high-Tc materials whose energy gaps vary by a factor of 30 . 

To motivate the analysis below that firmly establishes the mode interpretation 
of of the peak/dip/hump spectra and its connection with neutron data, we note 
that the spectral dip represents a pairing induced gap in the incoherent part of 
the spectral function at (tt, 0) occurring at an energy A -\- Qq, where A is the 
superconducting gap and Qa is the mode energy. We can estimate the mode energy 
from ARPES data from the energy difference between the dip {A + Qq) and the 
quasiparticle peak {A). 
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Fig. 52. Doping dependence of the mode energy: (a) Spectra at (7r,0) showing the 
decrease in the energy separation of the peak and dip with underdoping. Peak and 
dip locations were obtained by independent polynomial fits and carefully checked 
for the effects of energy resolution, (b) Doping dependence of the collective mode 
energy inferred from ARPES together with that inferred from neutron data [tl2q]. 



In Fig. |52| b we plot the mode energy as estimated from ARPES for various 
doping levels as a function of Tc and compare it with neutron measurements. We 
find striking agreement both in terms of the energy scale and its doping depen- 
dence 1 125 1 . Th e sam e agreement, in greater detail, has been recently found using 
tunneling data jl2(j ], as shown in Fig. We note that the mode energy inferred 
from ARPES decreases with underdoping, just like the neutron data, unlike the 
gap energy (Fig, ^ofa), which increases. This can be seen directly in the raw data, 
shown in Fig. B2k. This is also seen from the tunneling data, where they have found 
that the mode energy scales with doping as STc, just like the neutron resonance. 
An interesting point from the tunneling is that the ratio of the mode energy to the 
gap energy saturates to 2 in the overdoped limit, as would be expected for a coUec- 
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tive mode sitting below a continuum with a gap of 2A. Moreover, there is strong 
correlation between the temperature dependences in the ARPES and neutron data. 
While neutrons see a sharp mode only below Tc, a smeared out remnant persists 
up to T* 127 1 . As the sharpness of the mode is responsible for the sharp spectral 
dip, one then sees the correlation with ARPES where the dip disappears above Tc, 
but with a remnant of the hump persisting to T* . 

An important feature of the neutron data is that the mode only exists in a 
narrow momentum range about {tt,-k), and is magnetic in origin. To see a further 
connection with ARPES, we return to the results of Fig. ^ Note the dispersion 
along the two orthogonal directions are similar (Fig. |49| ;) , unlike the dispersion in- 
ferred in the normal state. As these two directions are related by a (tt, tt) translation 
((a:,0) = {0, —x); {0, —x) + (7r,7r) — (7r,7r — x)), we see that the hump dispersion 
is clearly reflecting the (tt, tt) nature of the collective mode. This dispersion is also 
consistent with a number of models in the literature which identify the hi gh en ergy 
feature in the pseudogap regime as a remnant of the insulating magnet |128|. We 
note, though, that the mode is due to quasiparticle pair creation and thus not just 
a continuation of the spin wave mode from the antiferromagnet |129|. 

This brings up a question that is at the heart of the high Tc problem: how can 
a feature which can be understood as a strong coupling effect of superconductivity, 
as discussed above, turn out to have a dispersion that resembles that of a magnetic 
insulator? The reason is that the collective mode has the same wavevector, {n,n), 
which characterizes the magnetic order of the insulator. It is easy to demonstrate 
that in the limit that the mode energy goes to zero (long range order) , one actually 
reproduces a symmetric dispersion similar to that in Fig. p9|c, with the spectral gap 
determined by the strength of the mode ||4l|. This is in accord with the increase 
in the hump energy with underdoping (Fig. |5C|b) tracking the rise in the neutron 
mode intensity. Since the hump scales with the superconducting gap, the obvious 
implication is that the mode is intimately connected with pairing, a conclusion 
which c an a lso be made by relating the mode to the superconducting condensation 
energy [L30]. That is, high Tc superconductivity is likely due to the same magnetic 
correlations which characterize the insulator and give rise to the mode. 



7.7 Dispersion Kink of Nodal Quasiparticles 

So far, our discussion has largely been centered on behavior near the (7r,0) point 
of the zone. We now turn to consideration of other k vectors. 

Remarkably, we find that the effects discussed above are manifest even on the 
zone diagonal where the gap vanishes, with significant changes in both the spectral 
lineshape and dispersion below Tc, relative to the normal state where the nodal 
points exhibit quantum critical scaling ji^]. Specifically, below Tc a kink in the 
dispersion develops along the diagonal at a finite energy (~70 meV) |^,^. This 
is accompanied, as required by Kramers-Kroni g re lations, by a reduction in the 
linewidth leading to well-defined quasiparticles P3|. As one moves away from the 
node, the renormalization increases, and the kink in dispersion along the diagonal 
smoothly evolves into the spectral dip, with the same characteristic energy scale 
throughout the zone. 

In Fig. we plot the dispersion of the spectral peak above Tc obtained from 
constant k scans (energy distribution curves or EDCs), and the peak in momentum 
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Fig. 53. 
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obtained from constant uj scans (momentum distribution curves or MDCs) from 
data for a rc=90K sample along the (7r,7r) direction |^^. We find that the EDC 
and MDC peak dispersions are very different, a consequence of the uj dependence 
of E. 

To understand this, we start by noting that since ~ v%{k — kp), then from 
Eq. 8 the MDC at fixed w is a Lorentzian centered at k — kjr + [uj — X!'{lo)] /v%, 
with a width (HWHM) Wm = \S" {lo)\/v%, provided (i) S is essentially independent 
of k normal to the Fermi surface, and (ii) the dipole matrix elements do not vary 
significantly with k over the range of interest. That these two conditions are fulfilled 
can be seen by the nearly Lorentzian MDC lineshape observed in ARPES 

On the other hand, in general, the EDC at fixed k has a non-Lorentzian line- 
shape refiecting the non-trivial oj-dependence of E, in addition to the Fermi cutoff 
at low energies. Thus the EDC peak is not given \yy u) ~ v%{k — kp) — E'{u)) — 
but also involves E" , unlike the MDC peak. Further, if the EDC peak is sharp 
enough, making a Taylor expansion we find that its width (HWHM) is given by 
We ^ \E"{Ek)\/[l - dE' /du\E^], where Ek is the peak position. 

We see that it is much simpler to interpret the MDC peak positions, and thus 
focus on the change in the MDC dispersion going from the normal (N) to the 
superconducting (SC) state shown in Fig. The striking feature of Fig. Isij a is 
the development of a kink in the dispersion below Tc. At fixed uj let the dispersion 
change from fcjv to ksc- Using v%{kN — ksc) = ~ '^iv('^)i directly obtain 

the change in real part of E plotted in Fig. B4c. The Kramers-Kronig transformation 
of E'sc - E'j^ then yields - Egc, plotted in Fig. which shows that \Esc\ 
is smaller than \E']^\ at low energies. 

We compare these results in Fig. pBj a with the Wm = \ E"\/vp estimated directly 
from the MDC Lorentzian linewidths. The normal state curve was obtained from 
a linear fit to the corresponding MDC width data points in Fig. |55|a, and then 
the data from Fig. |54| d was added to it to generate the low temperature curve. 
We are thus able to make a quantitative connection between the appearance of a 
kink in the (MDC) dispersion below Tc and a drop in the low energy scattering 
rate in the superconducting state relative to the normal state, which leads to the 
appearance of quasiparticles below Tc [Q. We emphasize that we have estimated 
these T-dependent changes in the complex self-energy without making fits to the 
EDC lineshape, thus avoiding the problem of modeling the ui dependence of E and 
the extrinsic background. 

In Fig. ^sjb, we plot the EDC width obtained as explained in Ref. from 
Fig. ^sjd. As an interesting exercise, we present in Fig. ^sj: the ratio of this EDC 
width to the MDC width of Fig. [s^ a (dotted lines), and compare it to the renormal- 
ized MDC velocities, 1/v = dk/diu, obtained directly by numerical differentiation 
of Fig. (solid lines). We note that only for a sufficiently narrow EDC lineshape 
is the ratio We/Wm — Vp/[1 — dE' /dw] = vp. Interestingly, only in the supercon- 
ducting state below the kink energy do these two quantities agree, which implies 
that only in this case does one have a Fermi liquid. 

Similar kinks in the dispersion have been seen by ARPES in normal metals 
due to the electron-phonon interaction [131]. Phonons cannot be the cause here, 
since the kink disappears above Tc- Rather, this effect is suggestive of coupling to 
an electronic collective excitation which only appears below Tc. Recently, this view 
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has been challenged by Lanzara et a/. |132{ , and we discuss this work at the end of 
this subsection. 
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Fig. 55. (a) Comparison of change in E" obtained directly from the MDC widths 
(HWHM) to the one obtained from the dispersion in Fig. Md by using the Kramers- 
Kronig transform, (b) HWHM width obtained from EDCs shown in (d). Lines 
marked by fit are linear in normal state and linear/cubic in superconducting state. 
The data in (b) fall below the fits at low energies because of the Fermi cut-off of 
the EDCs. (c) Renormalized MDC velocity obtained from differentiating Fig. 
(solid lines), compared to the ratio We/Wm from (a) and (b). (e) Ratio of EDC 
dispersion slopes above and below the kink energy at various points along the Fermi 
surface (from middle panels of Fig. 



We now study how the lineshape and dispersion evolve as we move along the 
Fermi surface. An analysis similar to the above is possible, but more complicated 
due to the presence of an energy gap |^]. We will thus confine ourselves here to 
a general description of the data. In Fig. ^ we plot raw intensities for a series 
of cuts parallel to the MY direction (normal state in left panels, superconducting 
state in middle panels). We start from the bottom row that corresponds to a cut 
close to the node and reveals the same kink described above. As we move towards 
(tt, 0), the dispersion kink (middle panels) becomes more pronounced and at around 
kx—0.55 develops into a break separating the faster dispersing high energy part of 
the spectrum from the slower dispersing low energy part. This break leads to the 
appearance of two features in the EDCs, shown in the right panels of Fig. M. Fur- 
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ther towards (7r,0), the low energy feature, the quasiparticle peak, becomes almost 
dispersionless. At the (vr, 0) point, this break effect becomes the most pronounced, 
giving rise to the peak/dip/hump in the EDC. We note that there is a continuous 
evolution in the zone from kink to break, and these features all occur at exactly 
the same energy. 

The above evolution is suggestive of the self-energy becoming stronger as the 
(tt, 0) point is approached. This can be quantified from the observed change in the 
dispersion. In Fig. [sS^ we plot the ratio of the EDC dispersion slope above and 
below the kink energy at various points along the Fermi surface obtained from 
middle panels of Fig. |5^. Near the node, this ratio is around 2, but becomes large 
near the (tt, 0) point because of the nearl y dispersionless quasiparticle peak. A 
different behavior was inferred in Ref. |132|, but in their case, the cuts near (tt, 0) 
were perpendicular to ours, and thus not normal to the Fermi surface. 

The lineshape also indicates that the self-energy is larger near (7r,0), as is evi- 
dent in Fig. ^ Along the diagonal, there is a gentle reduction in S" at low energies, 
as shown in Fig. [ssj a and b, with an onset at the dispersion kink energy scale. In 
contrast, near the (tt, 0) point there must be a very rapid change in E" in order to 
produce a spectral dip, as discussed above. Despite these differences, it is important 
to note that these changes take place throughout the zone at the same characteristic 
energy scale (vertical line in Fig. ^7| ). 

As also discussed above, the (tt, 0) ARPES spectra can be naturally explained 
in terms of the interaction of the electron with a collective mode of electronic 
origin which only exists below Tc. It was further speculated that this mode was the 
neutron resonance. Here we have shown that dispersion and lineshape anomalies 
have a continuous evolution throughout the zone and are characterized by a single 
energy scale. This leads us to suggest that the same electron-mode interaction 
determines the superconducting lineshape and dispersion at all points in the zone, 
including the nodal direction. In essence, there is a suppression of the low energy 
scattering rate below the finite energy of the mode. Of course, since the neutron 
mode is characterized by a (tt, tt) wavevector, one would expect its effect on the 
lineshape to be much stronger at points in the zone which are spanned by (7r,7r), 
as observed here. 

A similar conclusion has been reached by Johnson et al. |l33|] , where they find 
that the kink energy scales with doping like the neutron resonance (Fig. [sS^), 
and that the temperature dependence of E' tracks that of the resonance intensity 
(Fig. Is^b). Moreover, they find that S' increases with underdoping (Fig. ^sja), much 
like that extracted from the peak/dip/hump lineshape at (7r,0). 

Detailed calculations which take into account the momentum dependence of 
the n eutron resonance give an extremely good description of the experimental data 
134 1 . These results can be understood by studying the Feynman diagram of Fig. 
The key point is that the neutron resonance has a finite width in momentum space, 
corresponding to a short correlation length of order 2 lattice constants. Because of 
this, there is now an internal sum over momentum in the diagram, which will be 
dominated by the flat regions of the fermionic dispersion around the (tt, 0) points. 
This means that structure in the electron self-energy will occur at an energy of 
i^max + ^res, independent of external momentum. This explains why the energy 
scale is invariant throughout the zone. On the other hand, as the external momen- 
tum is swept, the momentum dependence of the neutron form factor is probed. 
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Fig. 58. Left panel: Plot of cuo, the energy of the maximum value of ReX" in the 
superconducting state (open circles), and ljq^' (gray circles), the energy of the max- 
imum in difference between the superconducting and normal state values plotted as 
a function of Tc referenced to the maximum (= 91 K). The coupling constant 

A (black triangles) is referenced to the right-hand scale. Right panel: Temperature 
dependence of KeE{LjQ^') from the nodal line for the UD69K sample (black squares) 
compared with the temperature dependence of the intensity of the resonance mode 
observed in inelastic neutron scattering studies of underdoped YBa2Cu306-i-a;, Tc 
= 74 K (gray circles) (adapted from Ref.|133|). 



Since the latter peaks at Q = (7r,7r), then the magnitude of the self-energy will 
be maximal at (7r,0), since these points are connected by Q, and minimum at the 
node. This explains why the peak/dip/hump effect first weakens into a "break" 
effect and then into a dispersion kink as the node is approached. Thes e ca lculations 
have been recently extended to incorporate bilayer splitting effects |135|, and are 
able to explain a number of unusual lineshape and dispersion features present in 
data on heavily overdoped Bi2212 

As mentioned above, this picture has been challenged by Lanzara et a/. |l32| . 
These authors claim that the kink is still present above Tc, except it is smeared 
in energy. Moreover, they find that a kink is present in a large variety of cuprates. 



98 



J.C. Campuzano, M.R. Norman, and M. Randeria 



including Bi2201 and LSCO, with an energy which is material and doping indepen- 
dent, as shown in Fig. |9| They argue that all of these observations are in support 
of a phonon interpretation of the kink. 

Although initially attractive, there are some problems with this scena rio. First, 
in regards to the kink above Tc, it has been claimed by Johnson et a/. |l33| that 
the "kink" above is simply the curvature in the dispersion one expects based on 
marginal Fermi liquid theory. In support of this, they argue that the maximum in 
the real part of X' is at a different energy in the normal state than in the supercon- 
ducting state. Though this appears to be the case, there is indeed residual structure 
in the normal state self-energy at the kink energy in the optimal doped sample we 
have looked at. On the other hand, our normal state data actually corresponds to 
the pseudogap phase, and as a residual of the neutron resonance is present in the 
pseudogap phase, the residual "kink" above Tc (if really there) does not rule out a 
magnetic intepretation. Moreover, as Johnson et aZ. convincingly showed (Fig. |58|b), 
there is definitely a large component to the self-energy which follows the same or- 
der paramter like temperature variation that the neutron resonance intensity does. 
This observation is supported by recent work of Gromko et oZ. [136| concerning a 
dispersion kink in the bonding band of heavily overdoped Bi2212 near (vr, 0). 

Second, in regards to constancy of the energy scale, this is indeed an inter- 
esting observation, though we note this statement contradicts that of Johnson et 
oZ.concerning the doping dependence of the kink energy mentioned above (Fig. [s^ ). 
Also, even in a phonon model, the kink energy should occur at the sum of the 
maximum gap energy plus the phonon energy. Why the sum of these two numbers 
should be doping and material independent is a real puzzle (as it would be for a 
magnetic interpretation as well). And, why only one phonon would be relevant, 
despite the large number of phonons present in the cuprates, is another puzzle. 
Still, a ph onon model for the kink has certain attractions, as discussed by Lanzara 
et a/. 1 132 1 . Certainly, more work is needed to definitively resolve the controversies 
surrounding the origin of the dispersion kink. 



7.8 Condensation Energy 

We conclude this section by discussing the relation of ARPES data to the super- 
conducting condensation energy. 

We begin with the assumption that the condensation energy does not have a 
component due to phonons. To proceed, we assume an effective single-band Hamil- 
tonian which involves only two particle interactions. Then, simply exploiting stan- 
dard formulas for the internal energy U = {H — ^N) {fi is the chemical potential, 
and A'^ the number of particles) in terms of the one-particle Green's function, we 
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Un-Us = 

J2 du;(Lo + ek)f{io)[AM{K^)-As{'k,'^)] (21) 



k 



where the spin variable has been summed over. Here and below the subscript A'^ 
stands for the normal state, S for the superconducting state. A{k,uj) is the single- 
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Fig. 59. Ubiquity of a sudden change ('kink') in the dispersion. Top panels are plots 
of the dispersion (derived from the momentum distribution curves) along (0, 0) — 
(tt, tt) (except panel b inset, which is off this line) versus the rescaled momentum k' 
for different samples and at different doping levels, a-c, Doping (5) dependence of 
LSCO (at 20 K; a), Bi22f2 (superconducting state, 20 K; b), and Bi220f (normal 
state, 30 K; c). Dotted lines are guides to the eye. The kink position in a is compared 
with the phonon energy a.t q = (vr, 0) (thick red arrow) and the phonon width and 
dispersion (shaded area) from neutron data. The doping was determined from the 
Tc versus doping universal curve. Inset in b, dispersions off the (0, 0) - (tt, tt) 
direction, showing also a sharpening of the kink on moving away from the nodal 
direction. The black arrows indicate the position of the kink in the dispersions. 
d,e, Temperature dependence of the dispersions for LSCO (d, optimally doped) 
and Bi22I2 (e, optimally doped), f, Doping dependence of A along the (0, 0) - 
(tt, tt) direction as a function of doping. Data are shown for LSCO (filled triangles) 
and NdLSCO (1/8 doping; filled diamonds), Bi220I (filled squares) and Bi2212 
(filled circles in the first Brillouin zone, and unfilled circles in the second zone). 
The different shadings represent data obtaine d in d ifferent experimental runs. Blue 
area is a guide to the eye. (adapted from Ref.|132l). 
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particle spectral function, f{uj) the Fermi function, and the bare energy disper- 
sion which defines the kinetic energy part of the Hamiltonian. Note that the ^A'^ 
term has been absorbed into uj and et , that is, these quantities are defined relative 
to the appropriate chemical potential, fiN or ^s- In general, and /is will be dif- 
ferent. This difference has to be taken into account, since the condensation energy 
is small. 

The condensation energy is defined by the zero temperature limit of Un — Us 
in the above expression. Note that this involves defining (or somehow extrapolating 
to) the normal state spectral function at T = 0. Such an extrapolation, which we 
return to below, is not specific to our approach, but required in all estimates of 
the condensation energy. We remark that Eq. |2l| yields the correct condensation 
energy, A''(0)Z\^/2, for the BCS theory of superconductivity. 

We also note that Eq. can also be broken up into two pieces to individually 
yield the thermal expectation value of the kinetic energy (using 2ek in the paren- 
theses in front of f{iL>)), and that of the potential energy (using ui — instead). 

The great advantage of Eq. ^ is that it involves just the occupied part of 
the single particle spectral function, which is measured by ARPES. Therefore, in 
principle, one should be able to derive the condensation energy from such data, 
if an appropriate extrapolation of the normal state spectral function to T=0 can 
be made. On the other hand, a disadvantage is that the bare energies, e^, are 
a priori unknown. Note that these are not directly obtained from the measured 
ARPES dispersion, which already includes many-body renormalizations. Rather, 
they could be determined by projecting the kinetic energy operator onto the single- 
band subspace. 

Some of the problems associated with an analysis based on experimental data 
can be appreciated. First, the condensation energy is obtained by subtracting two 
large numbers. Therefore, normalization of the data becomes a central concern. 
Problems in this regard when considerin g n( k), which is the zeroth moment of the 
ARPES data, were discussed previously |138|. For the first moment, these problems 
are further amplified due to the uj weighting in the integrand. When analyzing real 
data, we have found that the high energy tail contribution to the first moment is 
very sensitive to how the data are normalized. Different choices of normalization 
can even lead to changes in sign of the first moment. 

Another concern concerns the k sum in Eq. |2^. ARPES has k-dependent matrix 
elements, which lead to weighting factors not present in Eq. |2^. These effects can 
in pr inciple be factored out by either theoretical estimates of the matrix elements 
1 139 1, or by comparing data at different photon energies to obtain information on 
them §. 

Another issue in connection with experimental data is an appropriate extrapo- 
lation of the normal state to zero temperature. Information on this can be obtained 
by analyzing the temperature dependence of the normal state data, remembering 
that the Fermi function will cause a temperature dependence of the data which 
should be factored out before attempting the T — extrapolation. We finally note 
that the temperature dependence issue is strongly coupled to the normalization 
problem mentioned above. In ARPES, the absolute intensity can change due to 
temperature dependent changes in absorbed gasses, surface doping level, and sam- 
ple location. Changes of background emission with temperature is another potential 
problem. 
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Despite these concerns, we believe that with careful experimentation, many of 
these difficulties can be overcome, and even if an exact determination of Eq. 
is not possible, insights into the origin of the condensation energy will certainly 
be forthcoming from the data. This is particularly true for ARPES, which has 
the advantage of being k resolved and thus giving one information on the relative 
contribution of different k vectors to the condensation energy. 

Insights into what real data might indicate have been offered by us [137| in 
the context of the "mode" model illustrated in Fig. What we found was that 
for parameters characteristic of optimal doped ARPES data, the superconducting 
condensation was driven by kinetic energy lowering, as opposed to the potential 
energy lowering found in BCS theory. This occurs because Uk becomes sharper in 
the superconducting state than in the normal state. In essence, the normal state is 
a non Fermi liquid and the superconducting state is a Fermi liquid, so what occurs 
is that the effect of quasiparticle formation on sharpening rik is greater than the 
effect of particle-hole mixing on smearing it. The net result is a sharpening, leading 
to a lowering in kinetic energy. In BCS theory, the normal state is a Fermi liquid, 
and thus only the particle-hole mixing effect is present, leading to a net smearing 
of rife and thus an increase in the kinetic energy. The same model can be used to 
evaulate the optical sum rule [14C], and what is found is a violation of t he s um 
rule with a sign and magnitude consistent with recent optics experiments |141|. It 
will be of great interest to see whether these results can be confirmed directly from 
ARPES data, as speculated early on by Anderson [142]. 
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